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DIE ROLLE DER BEUGUNGSWELLE ІМ DEN 
FRAUNHOFERSCHEN BEUGUNGSERSCHEINUNGEN 


von A. RUBINOWICZ 


Staatliches Mathematisches Institut, Warszawa 
(Eingegangen am 10. April 1953) 


Die Fraunhoferschen Beugungserscheinungen lassen sich im Rahmen der Kirch- 
hoffschen Theorie der Beugung durch die Interferenz der von den einzelnen Elementen 
des beugenden Randes ausgehenden Beugungswellen erklären. Die „Fraunhofersche 
Beugungswelle" kann dabei durch einen einfachen Grenzübergang aus dem vom Ver- 
fasser (1917) für den Fresnelschen Fall angegebenen Ausdruck gewonnen werden. Dabei 
kann der beugende Rand auch durch eine Raumkurve gegeben sein. Im Gegensatz zu 
den Fresnelschen Beugungserscheinungen kommt im Falle einer ganz im Endlichen 
liegenden Beugungsöffnung die direkt einfallende Lichtwelle in der gesamten Beugungs- 
erscheinung gar nicht speziell zur Geltung. Der angegebene Ausdruck für die Fraunho- 
fersche Beugungswelle erleichtert die Behandlung von Spezialfällen dieser Beugungser- 
scheinungen, vornehmlich in dem Falle, wo der beugende Rand aus einer Raumkurve 
besteht. 


$1. Problemstellung 


Die Wellenfunktion, die in der Kirchhoffschen Theorie die Fresnelschen 
Beugungserscheinungen beschreibt, lässt sich in eine direkt einfallende und 
eine Beugungswelle aufspalten. Die direkt einfallende Welle wird in allen Raum- 
punkten, in denen nach der geometrischen Optik im Falle eines nichtreflektieren- 
den und undurchsichtigen Schirmes Licht vorhanden ist, durch die Funktion 
gegeben, die die unbehinderte Ausbreitung des Lichtes beschreibt. In allen übri- 
gen Raumpunkten aber verschwindet sie. Die Beugungswelle stellt eine Wellen- 
bewegung dar, die im Sinne der Youngschen Anschauungen aus den an den 
einzelnen Elementen des beugenden Randes erzeugten Sekundärwellen besteht. 

Im Falle der Fraunhoferschen Beugungserscheinungen rücken sowohl die 
Lichtquelle als auch der Beobachtungspunkt unendlich weit vom beugenden 
Schirme weg, so dass der Lichtkegel der einfallenden Welle in eine exakt bestimmte 

Richtung zusammenschrumpft. Man erwartet daher, dass die gesamte Fraun- 
hofersche Beugungserscheinung, mit eventueller Ausnahme der zur Einfalls- 
richtung des Lichtes entgegengesetzten Blickrichtung, durch Interferenz der 
verschiedenen Teile der Beugungswelle zustande kommt. 

Schwierigkeiten bereitet bei dem angegebenen Grenzübergang nur das Ver- 
halten der direkt einfallenden Lichtwelle, da ihre Richtungsverteilung vollkommen 
unstetig wird. Sie kann nur in der zur Einfallsrichtung des Lichtes entgegen- 
gesetzten Blickrichtung von Null verschieden sein und verschwindet für alle 


(5) 
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anderen Beobachtungsrichtungen. Es wird im folgenden zu erértern sein, ob 
diese Unstetigkeit durch eine entsprechende Unstetigkeit der Richtungsabhän- 
gigkeit der Beugungswelle kompensiert wird oder nicht. 

Mit der oben ausgesprochenen Anschauung, dass die Fraunhoferschen 
Beugungserscheinungen durch die Beugungswelle verursacht werden, steht der 
sehr elegante und einfache Nachweis von Laue (1936) im Einklang, dass das 
über die Beugungsöffnung erstreckte Fraunhofersche Integral mit Hilfe des 
zweidimensionalen Greenschen Satzes in ein längs des beugenden Randes ver- 
laufendes Linienintegral umgeformt werden kann. 

Wir wollen uns von der Richtigkeit der obigen Behauptung über die Rolle 
der Beugungswelle ganz direkt in der Weise überzeugen, dass wir in dem Aus- 
druck für die Kirchhoffsche Beugungswelle den bekannten Grenzübergang aus- 
führen, durch den aus dem Kirchhoffschen das Fraunhofersche Beugungsintegral 
gewonnen wird und sodann zeigen, dass der in dieser Weise erhaltene Ausdruck 
für die Fraunhofersche Beugungswelle die durch das Fraunhofersche Integral 
gegebene Wellenbewegung liefert. Allerdings werden wir diesen Nachweis nur 
in dem Falle erbringen, wo die Beugungsóffnung ganz im Endlichen liegt. 


82. Die Beugungswelle 


im Falle der Fraunhoferschen Beugungserscheinungen 


Um zunächst die Kirchhoffsche Beugungswelle ug in die Fraunhofersche 
überzuführen, legen wir sie unseren Betrachtungen in der Gestalt zugrunde, in 
der sie vom Verfasser (1917) angegeben wurde. Sie kann vektoriell in der Form 
1 f е%зӘ (rx о) ds 
ив = —— ~ 


4л s ro ro tro 


(1) 


geschrieben werden, Hier bedeutet ds ein Linienelement des beugenden Randes B 


und die Vektoren 0 und r r legen die Lage des Randelementes ds in in bezug auf die 


> 


Lichtquelle L bzw. den Beobachtungspunkt P fest, so dass о = Lds, r — Pds ist. 
Die Umlaufsrichtung bei der Integration über den beugenden Rand B ist dabei 


durch die Rechtsschraubenregel der der Lichtquelle zugewandten Normale n an 
die Fläche F zugeordnet, die die beugende Öffnung verschliesst. Diese Fläche 
soll im folgenden kurz auch als die beugende Üffnung bezeichnet werden. Bei 
der Ableitung des Ausdruckes (1) wird schliesslich vorausgesetzt, dass die von 
der Lichtquelle L "aic E Wellenbewegung bei der unbehinderten Ausbrei- 
tung des Lichtes durch е” /о* gegeben wird, wo 0% die о „Lichtquelle L— 
Beobachtungspunkt P“ bezeichnet. 

Um in (1) den Grenzübergang zu den Fraunhoferschen Beugungserschei- 
nungen durchzuführen bezeichnen wir mit O einen fixen, etwa in der Nahe des 
beugenden Randes B oder auf ihm selbst liegenden Punkt. Die Richtungen, in 
denen die Lichtquelle L und der Beobachtungspunkt P von O aus gesehen i ins 


Unendliche rücken, sollen durch die beiden Einheitsvektoren La bzw. eae 


Fraunhofersche Beugungswelle 


Ux 


gegeben sein. Der Einheitsvektor a г bezeichnet somit die Richtung in der das 
Licht aus dem Unendlichen auf den beugenden Schirm auffällt und der Einheits- 
vektor ap die Richtung in der wir auf den beugenden Schirm blicken. 

Seien ọọ bzw. го die Entfernungen der Lichtquelle L bzw. des Beobachtungs- 
punktes P von dem Punkte O, so kónnen wir in der in (1) auftretenden Exponen- 


tialfunktion 
055 0 Ка, , r=r + Rap (2) 


setzen, falls R den Vektor bezeichnet, der die Lage des Randelementes ds in 
bezug auf den Punkt O festlegt. Ausserhalb der Exponentialfunktion genügt 
es aber im Integranden der Beugungswelle (1) 

Ө —— оа = 10. Gp (3) 
anzunehmen. Um die Amplitude und Phase der Wellenbewegung beim Grenz- 
übergang im Beobachtungspunkte P festzuhalten, dividieren wir überdies die 
Beugungswelle (1) durch den Faktor 

ей (rot e) 


(4) 
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Wir erhalten so für die Beugungswelle u, im Falle der Fraunhoferschen Beugungs- 
erscheinungen unmittelbar den Ausdruck 


: 1 - x a It > > 4 

ив = Е ur СРЕ (apt epu. (5) 
4л 1 -- ар ay, $ 

der nur von der Einfallsrichtung des Lichtes a, und der Beobachtungsrichtung 


ap abhängt. Bei der Ableitung von (5) wurden über die Gestalt des beugenden 
Randes B keine einschränkenden Voraussetzungen gemacht. Er kann somit auch 
aus einer Raumkurve bestehen. 

Der Ausdruck (5) soll im folgenden als die Fraunhofersche Beugungswelie 
bezeichnet werden. In dem Falle, wo der beugende Rand B ganz im Endlichen 
liegt, ist er als Funktion der Blickrichtung аһ überall endlich und mit Ausnahme 
der Blickrichtung ap = ae , die zur Fortschreitungsrichtung des einfallenden 
Lichtes entgegengesetzt gerichtet ist, auch überall stetig. In der zuletzt genannten 
Richtung ist aber nur sein Grenz- und nicht auch sein Funktionswert definiert. Um 
uns davon zu überzeugen bezeichnen wir mit y den Winkel zwischen den beiden 


Vektoren ap und as Der Absolutbetrag des in dem Ausdruck (5) vor dem Inte- 
gral auftretenden Richtungsfaktors wird dann durch 
Ra к ылу ер. | (9 
1+ apa; 1 + cos g 2 | 
gegeben. Er wird also unendlich wenn wir nach der Lichtquelle blicken, d.h. 


wenn ар = —a, und daher p = л wird. 
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Dennoch bleibt in diesem Grenzfall die Fraunhofersche Beugungswelle 
endlich, da dann auch das in (5) auftretende Integral verschwindet. Mit Rück- 
sicht darauf, dass für jede geschlossene Kurve f ds =0 ist, kann nämlich das 
Integral in der Gl. (5) auch in der Form 


(еее С) а 1) ds (7) 


B 


dargestellt werden. Dieser Ausdruck verschwindet aber im Grenzfalle ар > — a, 


+ E - > Ф , 
proportional zu |ap + a,| = | / (ap + aj)? = 2 cos ^ Die Beugungswelle 


(5) ist somit im Falle der Fraunhoferschen Beugung in jeder Blickrichtung ap 


endlich und mit Ausnahme der Blickrichtung ap = —a, auch stetig. 

Durch die obigen Überlegungen wird jedoch nicht die Möglichkeit ausge- 
schlossen, dass der Ausdruck (5) in dem Grenzfall ар — I unendlich wird, wenn 
wir in ihm die Integration nur über einen Teil des beugenden Randes ausführen, 
weil dann das Integral nicht proportional zu | tp + "m verschwindet. Dass ein 
‚ solches Unendlichwerden tatsächlich eintritt, hat Laue (1936) für ebene Beugungs- 
ründer gezeigt. Y % 

Das Verhalten der Fraunhoferschen Beugungswelle (5) im Grenzfall ap > —a, 
wird durch die Tatsache verständlich, dass die Beugungswelle (1) im Falle der 
Fresnelschen Beugungserscheinungen in der Schattengrenze unstetig ist. 


$3. Vergleich der Fraunhoferschen Beugungswelle 


mit dem Fraunhoferschen Beugungsintegral 


In § 2 haben wir gesehen, dass die Fraunhofersche Beugungswelle (5) in der 
zur Fortschreitungsrichtung der einfallenden Welle entgegengesetzten Blickrich- 


tung ар = —a, einen endlichen Grenzwert aufweist. Um festzustellen, welche 
Rolle die einfallende Welle in den Fraunhoferschen Beugungserscheinungen 
spielt, wollen wir nun diesen Grenzwert mit dem Wert vergleichen, den das Fraun- 
hofersche Beugungsintegral für diese Blickrichtung liefert. Das letztere müssen 
wir dabei für den Fall nicht ebener beugender Rander zur Verfügung haben. 

Sei F eine Beugungsóffnung, die durch den beugenden Rand B begrenzt 
ist, so wird die Wellenfunktion in der Kirchhoffschen Theorie der Fresnelschen 
Beugungserscheinungen durch 


T e*t д ее | ei^? д ет 
шк Ye) Mr ee ee (8) 
: Е [ 


gegeben, wo п die der Lichtquelle L zugewandte Normale an die Flüche F bedeutet. 
Mit Rücksicht auf (2) und (3) geht dieser Ausdruck nach der Division durch 
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den Faktor (4) bei dem in $2 angegebenen Grenzübergang in das Fraunhofersche 
Beugungsintegral 


Шве x n (ag — ap) eltRGp + ар) df (9) 

über. 
Liegt die beugende Üffnung F ganz im Endlichen, so stellt (9) eine ausnahms- 
los für alle Blickrichtungen ap stetige Funktion dar. Für die spezielle Blickrich- 


tung ар = —a, ergibt das Fraunhofersche Beugungsintegral (9) insbesondere 


den Ausdruck . 
ا‎ пр 
Up аре a= 7220, К, (10) 


wo F den Vektor 


P= [sar (11) 
F 


bezeichnet. Aus der Tatsache, dass das Integral f n df für jede geschlossene 


Fläche verschwindet, folgt, dass der „Flächenvektor“ F, (11), unabhängig von 
dem übrigen Verlauf der Fläche F schon durch ihren Rand B eindeutig fest- 
gelegt wird. Seine Projektion auf die Normale einer Ebene gibt dabei den Flä- 
cheninhalt der Fläche an, den die Projektion der Kurve B auf dieser Ebene ein- 
schliesst. Der Umlaufssinn auf der Raumkurve B ist dem Richtungssinn des 
Vektors durch die Rechtsschraubenregel d.h. in der gleichen Weise zugeordnet, 


wie der Normalen n an die Fläche F. 
Um den Grenzwert der Fraunhoferschen Beugungswelle (5) für die Blick- 


richtung ap = —a, anzugeben, genügt es das Verhalten des Integrales (7) in 
Betracht zu ziehen, das dem in (5) auftretenden Integral gleich ist. Für kleine 


Werte von ap + ар wird es in erster Nüherung durch das Integral 
ik jJ [R( ap + а,)] ds (12) 
B 


gegeben. Um es zu berechnen benutzen wir ein System von rechtwinkligen, 
kartesischen Koordinaten Ё, 7, £, deren & -Achse die Richtung des Vektors ap + ор 


besitzt. Das Integral (12) wird dann gleich | 
ik(ap + aj) ¢ (J £ d£ + ға, + kf ED). (13) 


Das Integral f & d£, das über die &-Achse zu erstrecken ist, verschwindet. Die 


‘Integration in den Integralen f Е dn und f Е dC verläuft in der &,n bzw. £,£ -Ebene 


längs der Projektion В, bzw. В, дег Raumkurve, durch die der beugende Rand 


B gegeben wird. Unter der Annahme, dass die Projektion des Vektors F auf die 
£-bzw. n-Achse positiv ist, werden die Kurven В; und B, in dem in der Abb. 1 
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angegebenen Sinne durchlaufen. Da der Flicheninhalt, der von den Kurven В; 
und B, eingeschlossenen Flächen gleich F, bzw. F, ist, erhalten wir für das Inte- 
gral (13) schliesslich den Ausdruck 


il | [Ë (ap + ар] ds = ik (ap + a), G F, — EF) = 
B 
=: — ik (ap + aj) Х Ё. (14) 
Um zu егкеппеп,. dass das Fraunhofersche Beugungsintegral (5) im Grenz- 


falle ap > —a, einen Grenzwert annimmt, der gleich (10) ist, muss man noch 
bemerken, dass mit Riicksicht auf 


(ap + а)? = 2 (1 + ара), 


(ap + ar) X (aj — арр= 2 ap X а, 


die Fraunhofersche Beugungswelle (5) auch in der Gestalt 


ив = 


Ad a (ap + а) x ei ap) е Кар +21) ds (15) 
4m (ap + ар)? B 

dargestellt werden kann. Im Falle kleiner (ap + aj) -Werte kann für das hier auf- 
tretende Integral der Ausdruck (14) eingesetzt werden. Wir erhalten dann mit 
Rücksicht auf die Relation (a x b) (c = du (ac) (bd) — (ad) (bo) sowie darauf, 
dass (a pt aj) (оа ау = 0 ist, den Ausdruck (10). 


§4. Die Lauesche Umformung 


des Fraunhoferschen Beugungsintegrals 


Die Fraunhofersche Beugungswelle (5) muss man selbstverstündlich auch 
durch eine direkte Umformung des Fraunhoferschen Beugungsintegrals (9) 
erhalten. Es treten aber dann im Falle von beugenden Randern B, die durch 
Raumkurven gegeben werden, Schwierigkeiten von der gleichen Gróssenordnung 
auf, wie die, denen man bei der Herleitung des Ausdruckes für die Beugungs- 
welle (1) aus dem Kirchhoffschen Beugungsintegral (8) begegnet. 
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Wir wollen daher hier nur die Fraunhoferschen Beugungswelle (5) mit dem 
Ausdruck vergleichen, den Laue (1936) durch eine Umformung des Fraunhofer- 
schen Beugungsintegrals (9) in dem Falle gewonnen hat, wo der beugende Rand B 
durch eine ebene Kurve gegeben wird. Das in (9) auftretende skalare Produkt 


n(a, — ap) ist dann konstant und kann daher vor das Integralzeichen gesetzt wer- 
den. Um den Laueschen Ausdruck anzugeben, nehmen wir an, dass der Punkt O, 


in dem der Anfangspunkt des Vektors RI + jy liegt, sich in der gleichen 
Ebene E befindet, in der der beugende Rand B verläuft und bezeichnen mit h 


die als Vektor aufgefasste Projektion des Vektors ap + а; auf die Ebene Е. Es 
wird dann 


R(ap + a,) = Rh, | (16) 

so dass der Integrand in dem Ausdruck (9) in der Gestalt 

ф = gikR(ap + ар) — RR 
geschrieben werden kann. Da aber ф еше Lósung der zweidimensionalen Schwin- 
gungsgleichung 

9p م2‎ 
EA 

ox? ду 
ist, kann der Integrand ф in (9) durch — Ag/k?h? ersetzt und die Gleichung (9) 
mit Hilfe des Greenschen Satzes schliesslich in 


+ КА? = 0 


РЕ В (17) 


umgeformt werden. N bezeichnet dabei die in der Ebene E liegende áussere Nor- 
male an die ebene Randkurve B. Nimmt man an, dass ebenso wie in (5) die Um- 


laufsrichtung längs des beugenden Randes B der Normalen n an die Ebene E 
durch die Rechtsschraubenregel zugeordnet ist, so bilden die Vektoren N, ds 


und n ein Rechtssystem, so dass hx ds = п hy ds wird. Wir kónnen daher 
(17) vektoriell in der Gestalt 


мү ipee a tbe) [Oe (18) 
si Ал h2 ў 


darstellen. бе TU. 

Da der Vektor A als Projektion des Vektors ap + a; zugleich mit ihm ver- 
schwindet, so wird der in (18) vor dem Integral stehende Richtungsfaktor ebenso 
wie der entsprechende Faktor in (6) unendlich, falls die Blickrichtung entgegen- 
gesetzt gleich der Richtung des einfallenden Lichtes wird. Dennoch bleibt der 
ganze Lauesche Ausdruck (18) auch in diesem Grenzfalle für eine ganz im End- 
lichen gelegene Beugungsöffnung F endlich, wie bereits Laue klargestellt hat. 
Das ist auch daraus zu entnehmen, dass — wie im nächsten Paragraph bewiesen 


wird — der Ausdruck (18) nur einen Spezialfall des Ausdruckes (5) darstellt. 
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Die Tatsache, dass für die Blickrichtung ap = эмы пиг der Grenz- nicht 
aber auch der Funktionswert von ир durch (18) gegeben wird, ist insofern auf- 
fallend, als der Ausdruck (18) nur eine Umformung des Fraunhoferschen Beu- 
gungsintegrales (9) darstellt, das ausnahmslos für alle Blickrichtungen ap stetig 
ist. Diese Tatsache wird aber wohl verständlich, sobald man bedenkt, dass die 
Lauesche Umformung in dem in Rede stehenden Grenzfalle № >0 nicht mehr 
ausführbar ist, da ja dann der Integrand 9 — 1 wird und daher nicht mehr eine 
Lósung der Schwingungsgleichung darstellt. 


85. Vergleich des umgeformten Fraunhoferschen Integrals 
mit der Fraunhoferschen Beugungswelle 


Es steht zu erwarten, dass der Lauesche Ausdruck (18) im Falle eines ebenen 
beugenden Randes B mit der Fraunhoferschen Beugungswelle (5) identisch ist. 
Die beiden Ausdrücke sind in der Tat sehr áhnlich. Der für eine beliebige ráum- 
liche Randkurve B giltige Ausdruck (15) geht nämlich in den nur für eine ebene 
Randkurve verwendbaren Ausdruck (18) über, falls wir in (15) den Vektor аһ + а, 
formell durch seine Projektion auf die Ebene E d.h. durch den Vektor h ersetzen. 

Dennoch erfordert der Nachweis, dass die beiden Ausdrücke (15) und (18) 
für ebene Randkurven B einander gleich sind, einige Überlegung. Um ihn zu 
erbringen, führen wir mit Rücksicht auf (16) die Abkürzungen 


а = ар- ap, В = а,-- ар, 
> Eu ры а 5 19 
yafe kRh ds | ^59 ne ds ( ) 
B : 1 


ein. Die Bedingung für die Gleichheit von (15) und (18) lautet dann im Falle 
a0, h =0 


(6 x ya (Bx ph 
ae h? 

oder 
(B x у) (ah?— ha?) = 0. (20) 
Um das Bestehen der Beziehung (20) zu beweisen, bemerken wir zunächst, 
dass der Vektor ah? — hat senkrecht auf dem Vektor a steht. was man mit Riick- 
sicht auf ah = h? sofort feststellt, wenn man ihn skalar mit a multipliziert. Daher 
verschwindet das skalare Produkt dieses Vektors mit dem Vektor B X y, da ja 
dieser Vektor parallel ist zu a. Diese letztere Tatsache ergibt sich daraus, dass 


sowohl В als auch у senkrecht sind zum Vektor а. Dass а: wem 0 erkennt man 
unmittelbar aus der durch | (19) gegebenen Definition dieser beiden Vektoren. 


Um noch zu zeigen, dass ay — 0 ist, setzen wir a = h g a, 1» WO a, die senkrecht 


auf die Ebene E stehende Vektorkomponente von a bezeichnet. Da y in der 
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Ebene E liegt, ist a Ё ve 0, so dass man nur noch hy = 0 zu beweisen hat, um 
die Richtigkeit von ay — 0 zu bestätigen. Der Nachweis für hy = 0 ergibt sich 
aber unmittelbar, wenn man die Ebene 
E in parallele, zum Vektor h senkrecht 
stehende Streifen unterteilt. Wie aus 


der Abb. 2 ersichtlich ist, hat Rh und 
daher auch der Integrand exp i k Rh 
des Integrales (19), das den Vektor y 
definiért, für alle innerhalb eines sol- 
chen Streifens liegenden Randelemente 


ds den gleichen Wert. Da die Vektor- 
summe je zweier in einem solchen 
Streifen aufeinander folgenden Rand- 


elemente ds zum Vektor A senkrecht 


gerichtet ist, muss auch der Vektor y 
diese Richtung haben. 

Für einen ebenen Rand B ist damit 
der Nachweis für das Erfiilltsein der 
Beziehung (20) und somit auch für die 
Identität der beiden Ausdrücke (5) und Abb. 2 
(18) im Falle а #0, h +0 erbracht. 


8 6. Schlussbemerkungen 


Der durch (2) und (3) definierte Grenzübergang kann in der durch 
die Gl. (1) gegebenen Beugungswelle ug und in dem durch die Gl. (8) 
gegebenen Kirchhoffschen Integral nur dann einwandfrei ausgeführt werden, 
wenn der beugende Rand B bzw. die Beugungsóffnung F ganz im Endlichen 
liegen. Nur in diesem Falle werden die durch diesen Grenzübergang erhaltenen 
Integrale, nümlich die Fraunhofersche Beugungswelle (5) und das Fraunhofer- 
sche Beugungsintegral (9) im gewóhnlichen Sinne konvergent sein und nur in 
diesem Falle kann unser Resultat, dass die einfallende Welle sich im Beugungs- 
bild nicht direkt bemerkbar macht, als bewiesen angesehen werden. 

Sobald aber die Beugungsóffnung sich nur teilweise ins Unendliche erstreckt, 
geht auch der beugende Rand ins Unendliche und die Integrale in den Funktio- 
nen (5) und (9), die wir als die Fraunhofersche Beugungswelle bzw. ais das Fraun- 
hofersche Beugungsintegral bezeichnet haben, nehmen den Charakter von Dirac- 
schen 6-Funktionen an. Das im Fraunhoferschen Beugungsintegral ‚auftretende 
Integral wird direkt durch ein Produkt zweier ö-Funktionen 


+оо | 
ILE I өзек £35) dx dy = (22)! ò (kh,) ò (kh,) 
і 


gegeben, falls wir mit dem beugenden Rand allseitig ins Unendliche gehen. Da nun 
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zugleich mit dem Fraunhoferschen Beugungsintegral (9) auch die F raunhofersche 
Beugungswelle (5) sich wie eine Diracsche ó-Funktion verhilt, so besteht wiederum 
die Möglichkeit, dass die einfallende Welle sich nicht besonders bemerkbar macht. 

Dass die Wellenbewegung bei den Fraunhoferschen Beugungserscheinungen 


im Falle einer unendlich grossen Beugungsóffnung in der Blickrichtung ap— —a, 
wie eine 6-Funktion unendlich wird, ist verstándlich weil ja dann in dieser Rich- 
tung eine unendlich grosse Energie konzentriert wird. Die Frage, welches Aus- 
sehen, die in einem solchen Falle auftretende Beugungserscheinung insbeson- 


dere in der Nähe der Richtung m — a, besitzt (vgl. Boersch, 1952) kann jedoch, 
wie ich glaube, nur auf Grund einer Analyse der angewandten physikalischen 
Beobachtungsmethode entschieden werden. 

Wir bemerken aber, dass falls eine Fraunhofersche Beugungserscheinung 
wie gewóhnlich visuell durch ein Fernrohr beobachtet oder mit Hilfe eines Ob- 
jektivs auf einer photographischen Platte photographiert wird, die beugende 
Offnung und der beugende Rand stets ganz im Endlichen liegen, wie wir es bei 
den obigen Überlegungen vorausgesetzt haben. 

Man kann sich leicht an Beispielen überzeugen, wie gross die Vereinfachung 
ist, die bei der Verwendung der Fraunhoferschen Beugungswelle (5) zur prakti- 
schen Berechnung von Fraunhoferschen Beugungserscheinungen an beugenden 
Rändern erzielt wird, sobald diese durch Raumkurven gegeben werden. Das in (5) 
auftretende Integral kann z.B. sehr leicht für Gerade (vgl. Laue, 1936) berechnet 
werden, so dass für einen beugenden Rand, der durch einen ráumlichen Polygon- 
zug gegeben ist, die Fraunhofersche Beugungswelle durch eine Summe von ein- 
fachen Ausdrücken dargestellt wird. Da das Verhalten von Fraunhoferschen Beu- 
gungswellen, die von geradlinigen Begrenzungen ausgehen, leicht diskutiert wer- 
den kann, lässt sich auch unschwer das allgemeine Aussehen der in einem solchen 
Falle zu erwartenden Beugungsfigur voraussehen. 

Leider lässt sich die im Falle der Fresnelschen Beugungserscheinungen SOT 
nützliche Methode der stationären Phase (Rubinowicz, 1924) im Fraunhoferschen 
Fall nicht verwenden, da hier die „wirksamen Bereiche“ nicht von infinitesimaler 
Gróssenordnung sind. 

Bezüglich des in der Fraunhoferschen Beugungswelle (5) oder (15) auftre- 
tenden Integrals sei bemerkt, dass es einen Vektor darstellt, der senkrecht auf 


dem Vektor ap + a, steht. Einen Spezialfall dieser Behauptung stellt die im 
§5 bewiesene Tatsache dar, dass die Vektoren h und y, (19), senkrecht auf ein- 
ander stehen. Diese Behauptung kann ähnlich wie in 85 bewiesen werden, falls 
man stets zwei Elemente ds des beugenden Randes B zusammenfasst, die zwischen 


zwei parallelen und senkrecht auf dem Vektor a pt a, stehenden Ebenen liegen. 


КРАТКОЕ СОДЕРЖАНИЕ 


В. Рубинович, Роль дифракционной волны в теории дифракционных 
явлений Фраунгофера. 

В рамках теории Кирхгоффа имеем возможность “понять возникание 

дифракционных явлений Фраунгофера, предполагая что возникли они благо- 


EEE 
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даря интерференции дифракционных волн, распространяющихся OT частных 
элементов края дифракции. При этом, с помощью простого граничного пере- 
хода из формулы, предложенной автором в 1917 году для френельского 
случая, можем получить ,,Фраунгоферову дифракционную волну”. При этом, 
край дифракции может быть дан через пространственную кривую. Тогда 
как в дифракционных явлениях Френеля дифракционная волна играет 
важную роль, He обнаруживается она отчетливо B дифракционных явлениях 
Фраунгофера. Е 

Предложенная формула дифракционной волны Фраунгофера значительно 
облегчает дискуссию соответствующих, специальных примеров, прежде всего 
в случае, когда край дифракции дан через пространственную кривую. 
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ON A REGULAR FIELD THEORY II (QUANTIZED)* 


By Jerzy Rayskı 


Department of Physics, Nicholas Copernicus University, Torun 
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A non-local field theory constructed in a domain restricted by two hypersurfaces 
is quantized by a lagrangian method developed by Schwinger. The usual canonical commu- 
tation relations hold for pairs of points on either of the two hypersurfaces. No repre- 
sentations connected with intermediate surfaces exist. The Schrödinger equation does 
not exist. The two boundary hypersurfaces are physically distinguished by the fact that 
interactions with the apparatus of measurement take place on them. Probability ampli- 
tudes for transitions between states on these surfaces are computed. A method of a direct 
quantization of the field equations is discussed. Sufficient conditions for the existence 
of analytic solutions are given in the case of a quantized field coupled to an external field. 
The general case of two interacting quantized fields is discussed and some necessary 
conditions for the existence of regular solutions are formulated. 


1. Quasi-canonical quantization 


The non-local field theory derivable from the action functional 


w= | асе (х) ПЕС (1) 
2 Q 


in a space-time domain 0 restricted by two space-like hypersurfaces o} and оз, 
may be quantized by a lagrangian method developed by Schwinger (1951), pro- 
vided the lagrangian of interaction contains no derivatives of the field variables Ф". 
In this case we may introduce the conjugate momenta attached to the points on 


the surfaces оу and o, 


a (x) = n, nx), (2) 
where 7, is the unit vector normal to the surface at the point x,, and 
д М» д НЫ) 


а Nr eiumod: € a 3 
ла (x) = 29*(3 = 991) (3) 
We assume for simplicity that neither of the momenta vanishes identically, ven 
the system of fields is subject to no constraints. In consequence of the field equa- 
tions the variation of the action functional produced by varying the field quanti- 
ties depends only upon the values of др" on the boundaries с; and бу: 
LE و ب‎ 

* Part I appeared E this Journal Vol. XI, p. 314 (1953). Cf. also a paper in Il Nuovo 


i X, 1 (1953). 
Cimento ( ) (15) 
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à, W —(f— f) do a 0t. (4) 


(cf. Part I, formulae (1.5) and (3.2)). 


By taking a modified action functional 
py = Fili 9, (7, 9^) (5) 
we have the same equations of motion, so that W is equivalent to W. The variation 
of W is 
dW = àW — ô (f — f)do пар" = — (f—J)do дл" - 9* (6) 


and may be viewed as a variation of the action functionai produced by varying 
the л“’з regarded as original variables instead of the 045. This corresponds to 
a canonical transformation “> z^, л> — 0%. In this way some of the features 
of the canonical formalism appear also in the framework of this non-local theory. 

We assume that the general scheme of operators and state vectors in the 
Hilbert space holds also for a non-local field theory, and recapitulate briefly some 
of the formulae useful for a subsequent discussion. 

The infinitesimal unitary operator 


u-—1l-—iF, (7) 
where F is an infinitesimal hermitian operator, induces an alteration of the repre- 
sentation by changing the basic vectors into the new ones 

lau ыа аза (8) 
The representatives of an arbitrary vector | > and the matrix elements of an 
arbitrary operator G are altered by 


а) = аи >--<а| = Е], (9) 

|а иа —- <a >----і< ГІ ас>, (9°) 

ó«a'|G| a2 = <а |u1Gu|o" — «a'|G| a" 
zia | [Е.С] | a” >. (10) 


This infinitesimal change of the matrix elements produced by a change of the 
representation may be viewed alternatively as an (equivalent) alteration of the 
operator G itself (expressed in the unaltered representation) 

DEM Cad жаш (11) 


where 
| 6G — i[F, б]. (12) 
A simultaneous alteration of two basic systems of state vectors | a'>, Wits 
induced by infinitesimal unitary operators U, = l — iF, and U g = l — iF; pro- 
vides an alteration of the transformation function <a’ | В' >: 
à <a’ | "> = <а |u tu] B’>— <а | 8'> =i <a’ |F Е В". (13) 
We assume that in the non-local field theory the observables (i.e., physical 


quantities measured on o, or с») are hermitian operators constructed of field 
quantities taken at points lying either on the surface c, ог og. Thus, the state 
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vectors are | a’(c,)>, (or | B’(c.) >), where a’ (ог В’) denote n eigenvalu, of 
a complete set of commuting observables constructed on сү (or су). 

We notice that the variation Ô W given by (4) and the variation 6,W given 
by (6) are of the same form as the infinitesimal operator appearing to the right 
hand side of (13), so that a plausible assumption is that 6,W (and 6,W) describe 
infinitesimal alterations of the transformation function: 

д <a'(o) | B’(o.)> = iC < a'(o,) | ô W | B'(o) > 
= i C <a’ (o) | [do л° др — f do л° дф" | В (с) — (14) 


(and an analogous proportionality for д, W). ihe value of the dimensionless fac- 
tor C is a question beyond the scope of the formalism. To simplify the arguments, 
let us assume first that the field quantities are not spinors. In this case the in- 
tegro-differential field equations are of the second order in the derivatives, so 
that we may assume the values dy on c, and о, independently. By putting дф = 0 
on сз, the alteration (14) concerns the basic vectors on c, (and vice versa), whence 
the two terms in (14) separately play the róles of infinitesimal hermitian operators 
generating the changes of the representations on c, and o; respectively, 

Ел) = Cf do x“ др" , F,(oa) = C fdo п" ёр". (15’) 

о, Сз 

In the same way 


F,(o, = Cf do ê ng" , Е„(а) = С f do ôx" - p®. (15) 


As stated before, these alterations o the representations generated by {15) may 
be viewed equivalently as alterations of operato s. It is clear that these alterations 
are consequence of the infinitesimal changes дф“ (or дл“) of the field quantities. 
The change 

9" 9" Ч б" =u ‘и (16) 
is surely unitary if the variations дф" (or дл“) are c-numbers (i.e. commute with 
everything). We 1 mit ourselves to such variations only. 

The commutation relations between the field quantities y^ and л" taken at 
points lying on the same surface o, (or сз) follow from (12) and (15) with G — g* 
or G x 

òp (x) = iC | dr [° (x) p (9] 09^ ()), (17) 
01(0;) к 
ӧл? (x) = iC | do [л"(х), л? (x)] ôg" (x). (17°) 


Since 2° and p®, as nella orto and g* (x) at two different points x, x’ lying both 
on су (or both on 02) may be varied quite independently, we get 
l ГА 
[x*(x’), g^ (x)} = iC дав дао.) (-х jb | 
[7° (х), л? (x)] = 0. veumar(187) 
By taking F„ instead of F, we get in the same way 
; 19), v9] — 0, 


| 1 ! , " 
[v^(x)), л? (9] =— iC бав доз) (х — х ), (18 ) 
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where ЕЕ --ж) f (x') до = f(x) (19) 


for points x on c. Thus, the usual canonical commutation relations hold for both 
surfaces of measurements c, and oy if the factor C is unity (in natural units c 
= h = 1). Any other value of C would mean only a renormalization of the value 
of Planck's constant. 

In the case of spinor field quantities the equations are of the first order, but 
the field quantities p“ on the two surfaces су and оо, as well as the p%’s and z*'s 
on the same surface, may be varied as independent variables owing to the fact 
that the solution of a four-component Dirac equation involves four arbitrary 
functions, which corresponds precisely to the number of generating operators (15). 
Thus, we may apply also the above analysis to the case of spinor field quantities. 
As shown by Schwinger, the invariance under time reflection requires for spinor 
fields anticommuting variations др“ and дл“. This yields anticommutators in- 
stead of commutators in the formulae corresponding to (18). 

This procedure of quantization with the exclusion principle is self-consistent, 
since the infinitesimal transformations 9^ > 9°-+ дф“ with anticommuting 604 
are also unitary transformations, which may be seen from the fact that the new 
field quantities p*+ дф“ satisfy the same anticommutation relations as the old 
ones. 

Until now we have complete agreement between the local and the non-local 
quantum field theory, in particular the commutation relations (18) are identical 
with the well-known commutation relations provided by canonical quantization. 
Nevertheless it is hardly possible to call this method a "canonical quantization" 
on account of essential differences in comparison with the traditional theory. First 
of all, the commutation relations (18) apply to pairs of points lying on either of 
the boundary surfaces с; and с, but cannot be extended to pairs of points on any 
other space-like surface с situated between o, and og. Then, the energy and mo- - 
mentum operators P (c), and P, (o5) derived in Part I. do not satisfy the canonical 

a 
commutation relations [p*, Р] = — i a . The quantities Р (с) (where o is a 
p 
surface between c, and сз) are not constants of the motion and Р, (с) cannot be 
considered as a hamiltonian!. 

P,(c,) is not an observable in the usual meaning of this word, since it is not 
constructed of field quantities taken at points lying on the surface c, but (in con- 
sequence of the explicit appearance of the form-factor in P,(s,)) it involves the 
field quantities at all points of the domain Q. (Practically it involves only the 
points in a zone of the order of magnitude of the elementary time interval 4). 
Consequently a Schródinger equation does not exist in this formalism. This 
may be seen also from the fact that (in contradistinction to local field theory) the 


1 On account of the integral conservation laws P, (0 ) = P, (Ga), it is possible to find 
a quantity, attached to the surface g, which is constant and assumes the value P, (су) on с), 
but this quantity is not of the form of a surface integral. 


— Ө» 
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property of additiveness W, = Woo F Woo, 13 lacking owing to the threefold 


integrations appearing in (1). In the local theory it was possible to find a differential 
characterization of the time development by means of Wu. = Woo, + 0W 
with OW = Ро, but in the framework of the non-local formalism such a 


differential characterization is not possible. 

Inspite of the lack of differential characterization of the development of the 
system with time, we have an integral characterization by means of a unitary matrix 
connecting directly the results of measurements on the two surfaces су and о». 
Indeed, since the field quantities 0% and л“ satisfy the same commutation rela- 
tions on both surfaces c, and оз, the observables constructed by means of the same 
prescription on о; and 0; possess the same eigenvalue spectrum and are connected 
by means of a unitary transformation Иль: 


a(o,) = Uza (o) Из". (20) 
The transformation function connecting the two representations 
«a (o3) | a(o) > = <a’ (03) | Uiz’ | a” (оз) > (21) 


has the usual interpretation as probability amplitude. 

The connection between the matrices U,, and W,, may be seen in the follow- 
ing way. The fundamental assumption was that a change in the transformation 
function corresponding to a change of the field quantities by дф" is 


à <a’ (a) | a” (0) > = iC <a (o) | ôW | a (оз), (22) 


but from (21) the same change may be viewed also as an alteration of tlie unitary 


matrix (7/1; 
ô <a’ (o) | a” (с) > = <а (ea) | „Жир: | a^^ (o9) >. (23) 
An infinitesimal change of a unitary operator may be written as 
ОИ а, (24) 


where ôV is an infinitesimal hermitian operator. Hence 
ô <a'(o,) | a” (0) > = i <a (оз) | Ug Vie | а" (о) > 
=i <a’ (0) | Vy] «(e >. (25) 
From (22) and (25) it is seen that ОИ; is proportional to 6W,,. The connection 
between Л» and Wi, із 
Со, = U,,0U;. (26) 
The property of additiveness of the action functional was assumed hitherto 
on account of the composition law for transformation functions: 
< а, | а> ex scar | ag’ > dag" «ag |а" >, 
whence, by differentiation 
<a | SW yg | ast =f a | OW rs | as > dos" as" | as 
| Ја | as day" as" | | as 
or ôW, = ФР, + д зә compatible only with Wy» EU As + Wa, for any sur- 
face оз. The lack of additiveness in the non-local formalism shows that there exist 
no representations connected with intermediate surfaces оз. This is acceptable 


90 J. Rayski 


since the two surfaces o, and сз are ex definitione the only surfaces where meas- 
urements take place. Of course, we may perform measurements on оз instead of оз 
but, to describe this new situation, we have to use a new action functional WA, 
instead of Ито. 

The differences between the traditional local and this non-local field theory 
may be summarized as follows. In the local field theory the field equations are 
independent of the boundaries c, and сз which may be removed to plus and minus 
infinity. The field quantities are (formally) determined in the whole space-time mani- 
fold by the "initial" values on a given space-like surface c. This "initial" surface 
is not distinguished physically from any other space-like surface c'. The same 
solutions may be given with the aid of the "initial" values on o’ (coinciding оп со” 
with the solutions previously determined by given initial values on o). Thus, we 
may place the "initial" surface quite arbitrarily in the space-time (provided we 
assume suitable initial values on it) so that, in the traditional theory, the solutions ф° 
themselves acquire some sort of "absolute meaning", while the "initial" surface с 
is an auxiliary concept with a propositional character (since an interaction with 
external systems determining the initial state may occur on this surface or on 
any other surface). 

The situation in the non-local theory is different. The formalism is closely 
adapted to the fundamental physical problem: the initial state being prepared 
on сз, we ask for the state on another surface оу, where another interaction with 
the measuring apparatus actually takes place. The two surfaces are not abstractions 
and have not a propositional character, but are physically distinguished in conse- 
quence of disturbances provided by measurements. On the other hand, the field 
quantities and the equations of motion (field equations) acquire only a relative 
meaning with respect to the two surfaces of measurement. An immediate physical 
meaning is ascribed only to field quantities taken at points on the surfaces су and оу, 
while the o"'s at points lying between the two surfaces are without a direct relation 
to experiment and play only an auxiliary róle in the course of computations. 

Another significant difference between the local and this non-local formalism, 
is that the operators involving explicitely the form-factor, and therefore involving 
the field quantities at points lying between the surfaces of measurement (e.g., 
the energy of interaction) are not observables in the usual meaning of the word. 
This may, perhaps, facilitate a better understanding of some of the aspects within 
the theory of measurement, as e.g. the fact that we never measure directly the energy 
of interaction nor the total energy in a stationaty state, but infer its values indirectly 
from scattering experiments or experiments where transitions between two station- 
ary states are involved. 


ч 2. The matrix U, and the "direct quantization” 

The matrix U,, may be computed from the integral equations (2.8) of Part I 
as a power series of the coupling constant. We may write the field equations in an 
integral form using either о; ого, as the surface for which the initial values have 
been fixed 
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ф(х) = фз (x) + Í | j бахх" СФ, (x, и) FE 2 za (2.1) 
у 9g" (и), 


where с stands for o, ого, y^ denotes the hermitian conjugate of ф° (or the her- 


mitian conjugate multiplied by y, to the right in the case of fermions), u is one of the 
points x’,x’’,x’’’ which appear as argument of pf in L’, Фо, is the solution of the 
interaction-free wave equation (of Schródinger-Gordon or of Dirac) satisfying the 
same initial condition as p^ on су. If the surface o, is later than оз, G is identical 
with 44% (x — u) (ог 54%) and С} is identical with 4% (x — u) (or 574) for 
points x, u іп Q. Thus, the equations may be written as follows 


q^ (x) = фо (x) + ША Се) (x — и) hide (2:27) 
JJ 3g (u) 


[ es 
а(х) = 9? (x dx! dx! dx С (x — и) ——_, 2022 
вч = 95,0) + | [| 26-4 уу, 2.2") 

Q 
where С) stands either for 42” or 54), If we assume such initial values on 
бә which coincide with the solution of (2.2’) on o5, both solutions become identical 
and formulae (2.2) may be added or substracted from each other. If the initial 


function фо, is a general solution of the interaction-free equation, then ро is also 
a general solution of the interaction-free equation. On account of 


em + Git) NE р (24% _ Gre — G (2.3) 
we have | я 
9 P 
фо, (х) — Pa (x) = [| | f dx'dx''dx''' G(x — и) ЭРЫ) ; (2.4) 
24 IL? 
ф°(х) = ф(х) «fff dx'dx''dx'"' С (x — и) ap") ; (2.5) 
where 
1 
= +91) (2.6) 


С and С stand either for A and A or S and 5. | 
As shown in the foregoing section, the field quantities satisfy the usual canon- 
ical commutation relations for pairs of points on the surface о; (or сз). Therefore 
the general solutions of the interaction-free equations coinciding (together with 
their time-like derivatives in the case of bosons) with g^ on o, or су must satisfy 
the usual commutation relation for free fields, e.g. 
[v2 (а), pE] = Гр (а), v5, (| 
ud (x — x) if p’ is the complex conjugate of 9^ (27) 
РЗ 0 otherwise 
(corresponding anticommutation relations hold for fermions). The important P 
is that both sets of fields satisfy the same commutation relations, so that there 
exists a unitary operator U}, transforming Фо, into Фа: 
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pt = Шыр Um. (2.8) 
Гу is identical with the operator ( ` (20) of Section I. * rom (2.5) and (2.8), we have 


[U, pal = “a G(x — д 5 zo (2.9) - 
[U, фо] = „нае и) 5 : (2.9) 


The matrix U], may be computed as a power series in the coupling constant g 
either form (2.2)' together with (2.9’) or from (2.2'") together with (2.9”) or, 
most conveniently, from the symmetrized form of the field equations (2.4) together 
with half the sum of both equations (2.9) 


peur aL" 
Ноос 2 | г.|| ee — u) дода) | 4 (2.10) 


The method of solving the field equations (2.2) ог (2.4) (and formula (2.9) 
or (2.10)) by means of a power series expansion is a mechanical (automatic) pro- 
cedure and this fact ensures that the solutions will differ from those of the local 


theory by exactly the same substitutions which convert the corresponding formulae 
of local theory into the non-local ones. We may guess immediately the form of the 
operator U,,. It differs from the local one merely by the fact that every variable x 
appearing three times (either as an argument of the field quantities or as one of the 
arguments of the A-like functions) is replaced by three independent arguments 
x’, x", ж”, the corresponding term is multiplied by the form-factor F(x’, ж”, x’’’) 
and the integration f dx is replaced by threefold integrations / / f dx did dx ic 
This rule for an immediate transcription of the solutions (and of the operator U4.) 
into the non-local form is confirmed by the detailed computations in the appendix. 
The above described automatic substitutions cannot change the unitary character 
of the operator Uj, (the expression for U, obtained formally in the local theory 
is unitary) provided we use a form-factor guaranteeing the hermitian character 
of the lagrangian 3 
Fila Е а (2.11) 


This last argument was used in the previous papers of the author. The so called 
"direct quantization” introduced in those papers consists in assuming the usual 
commutation relations for the free waves (e.g. for ро). Then фе may be computed 
from the field equations and an operator U,, defined by Upo, = Po, Оа. Since 
this operator is a mechanical transcription of the analogous, well-known operator 
of the local theory, it must be unitary, whence Po, satisfies as well the usual commu- 
. tation relations for free waves. 

The arguments presented in Section 1 of this paper confirm the correctness 
of the "direct quantization" method in the case ise when the lagrangian of interaction 
contains no derivatives of the field quantities. It seems, however, that the method 
of "direct quantization" possesses a wider range of applicability and may be used 
К ани even in the cases when L contains derivatives and the ° "conjugated 
momenta” are essentially non-local: 
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a IK ӘР! 
л = ==; «ff с (2.12) 
9, де! 


For quantization it seems to be essential that the field equations are derivable from 
a variational principle. 


3. The problem of convergence 
We shall discuss the convergence problem in the non-local field theory on 
the basis of a simple model: a spinor field v coupled to a pseudoscalar field 9 
The interaction part of the action integral is 


а Ш dx! da! "du" P(x! 2", s") (рб), уэ vi) oi), GD 
Q 


where 


[a;b] =ab—ba , {a,b} = ab + ba. 


In order to simpliry the notations we shall write (3.1) symbolically by disregarding 
the proper order of the factors 


Же = ШЕУ x, и) (3.1”) 
9 


with the matrix y included in the definition of F. | 
Let us investigate first the comparatively simple case of а given external field ф. 
The field equations become linear integral equations for y and 9: 


ya) = ve(x) + 8 | dy N(x, y) p(y); (3,2) 
Q 
where w(x) is the interaction-free wave coinciding with the solution y(x) on ап 
initial surface (e.g. оз). The kernel N(x, y) is 


N(x, y) = || с (x, x") F(a’, x", y)o (x) dx’ ах”. (3.3) 
Q 
By a suitable choice of the form-factor the kernel is limited 


| N(y)|-«N or [venter (3.4) 
Q 


even in thecase of a singular external field ф with singularities of the type s How- 


ever, the value M depends upon the strength of the sources of the external field 
and for strong sources may be very large. The integral field equation (3,2) may 
be solved by iteration i 


v,(2) = vale) +8 7 dy N (s) v 6) (5) 


if the integral 4 


А- dy N(x, y) vol”) (3.6 
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exists, and if 
в || N(x, у) | dy < gM < 1 (we assume g > 0). (3.7) 


From (3.5) and (3.6) we have 


таса Sef rv n EI | vith [yol <E lAl 89 
whence 

| var 9, | < (gM)"g | 4| +> 0 (3.9) 

on account of (3.7), and the convergence of the sequence y to a solution follows. 


This solution may be written with the help of the resolving kernel R(x, y; g) as 
follows 


Ее | ITA NO (3.10) 
Q 
where 
REN CN ERE em) 
with 


М (asy) = fe [a N (xoa) Мена Ne: (3,12) 
2 2 
The series (3.11) converges if (3.7) is satisfied. (3.6) is a necessary condition, while 
(3.6) together with (3.7) are sufficient conditions for the existence of a solution 
expressible by means of a power series expansion in the coupling constant g. 
In a quantized theory equation (3.2) isasymbolic equation between operators, 
and may be replaced (in a representation) by an infinite set of equations for the 
matrix elements 


<a’ | p(x) | a^ = <a’ |ро(х) | a" > + ef N(x, y) <a’ |р (у) | a^ > (3.13) 
Conditions (3.6) becomes: the integrals 


<a’|A|a”> Jj dy N(x, y) <a’ | poly) | a” > (3.14) 

O. 
should exist for any states | a'>, | а'' >. In particular we may put for | а> 
a vacuum state | 0> and for <a’ | the state vector <0 | y(z) which yields a function 


—15-(у — z). Similarly one obtains a 15+(у — z) function. Thus, it is necessary 
that the following two integrals 


At = / dy N(x, y). S'o—2) (3.15) 


exist. As is well meus in the local field theory such integrals are divergent and 
give rise to infinite vacuum polarization effects. On the other hand, by a suitable 
choice of the form-factor, the integrals of type (3.15) may be made convergent. 
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1 


(Besides the removal of singularities of the type = the external field ф this is 
r 


just the chief reason for the introduction of a form-factor). The existence of the 
integrals (3.14) is guaranteed by the existence of (3.15), since an arbitrary solution 
of the interaction-free Dirac equation may be built up with the aid of the two funda- 
mental solutions 5+ and S- (ог S and S“). Thus, the existence of А+ is decisive 
for the case of a quantized field y interacting with an external field 9. On the other 
hand, condition (3. 7) (restricting the coupling constant, or the strength of the 
sources of the external field) may be abandoned and replaced by a less stringent 
one (3.4). As is well-known from the theory of linear integral equations, the resolv- 
ing kernel exists also in the case of an arbitrary value of the coupling constant g 
but, in general, is a quotient of two integer functions of the coupling constant. 
The solution is unique unless g assumes a value for which the resolving kernel 
possesses a pole. The existence of poles is surely possible if the initial surface og 
is removed to infinity. In this case the homogeneous equation obtained from (3.2) 
by dropping the free wave yg (x) is soluble and such solutions represent stationary 
states. On the other hand, the question whether (with a special choice of the ex- 
ternal field y) the resolving kernel R(x, y; g) possesses poles in the case of a proper 
domain Q restricted by two hypersurfaces is open?. 

At any rate, the non-local field theory with a suitable form-factor may be 
called a "regular field theory" since (for not too strong sources of the external 
field) the kernel is limited, the integrals (3.15) exist, and for sufficiently small 
values of the coupling parameter the solution is an integer function of g. 

The convergence problem in the case of two interacting fields (p not external) 
is much more involved. In order to have meaningful iteration formulae 


vH) = Yol) t Í | К, (х,у, 2) Pn (y) v. (2) dy dz (3.16) 
Q 


OR eae [ [| КИБИК: (8167) 
Q 


the following integrals must exist 


а= || акд. B= [ ®\] een Эм | 
Q Q 2 


(3.17) 
C= f dz | [ду К (у, фб), D=| I dy dz K (хуг) Фо(у)9 (2) 
Q Q Q 


and similar integrals E, F, G, H with the kernel K,. 
INNEREN 


3 The corresponding equation in local field theory is of a Volterra type so that poles seem to 
nt as (i) in the limit of a local theory the kernel N(x, у) 


be excluded. However, this is a weak argume өлі! | 
becomes singular (ii) the classification of the integral equations into those of a Volterra and of a Fred 


holm types looses its proper sense in the case of a non-local interaction. 
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In the appendix it is shown that, if the integrals A, B, ... H are all limited, 
say 
A.< М, BM THM (3.18) 
then the sequences y,, 9, converge to the solutions у, y for sufficiently small 
values of the coupling constant ku 
nen (3.19) 
4M 
In a classical theory the limited integrals (3.17) may be easily obtained by 
restricting the initial values of the field quantities, i.e. by restricting the free waves 


lg|< 


Yor Фо. On the other hand, in a quantum theory the yy sand фз are operators, and 
(in a representation) the matrix elements <a’ | yọ | a” > and <a’ | po | a^^ 
do not form a set of simultaneously restricted functions. Thus, the convergence 
proof given in the appendix cannot be applied directly to the case of quantized 
fields?. 

The classical restrictions upon the initial values of the field variables leading 
to inequalities (18) have a counterpart in the quantum theory in form of restric- 
tions upon the initial state. The author hopes that it will be possible to show that 
the iteration procedure converges if applied to expectation values of observables 
in states belonging to a certain class. The states of this restricted class will be 
denoted by | a’ >. The iteration procedure should also converge, if applied to 
computations of probability amplitudes for transitions between a state belonging 
to this class and any other state <a’ | U | a'> or <a” | U| a'>, where | а >> 
belongs to this class while | a’’> is an arbitrary state vector. The amplitudes 
<а | Ц |а’> should be integer functions of the coupling constant g, and the sum 


| [<a |U|a”’>da’<a” | U* | a^ (3.20) 
should be unity. Most probably this class of state vectors | a’ > may be character- 


ized by a limited expectation value of the kinetic energy E, of the system of fields 
on the initial surface: if a frame of reference exists in which 

<а | Е, |а’> <M = (3.21) 
then | a'> belongs to the class | a’>. The above supposition seems plausible 
on reasons of continuity (applicable to-the cases of integral field equations in 
a finite domain and with regular kernels). 

On the other hand, it seems improbable that the iteration procedure will 
yield satisfactory results in the general case where the initial state is quite arbi- 
trary. The following argument provides a hint that a field theory (even a regularized 
one) should possess extra solutions not expressible by means of the procedure of iter- 
CE uia c MUT Й ee op ER EE ee VIE MUST S UU 

3 Similarly to the case of a quantized ty coupled to an external @, the existence of the integrals 
B, C, D, F, G, H will be guaranteed for general matrix elements by the existence of similar integrals 
with St and /4= instead of Yo and фо. The integrals of the type B, С, F, С (with one function 5+ or 
A+) may be made limited in a finite domain by a suitable choice of the form-factor. The integrals 
of type D and Н (with two 5+ or one S+ and one 4+) are typical for the vacuum states. Such integrals 


may be compensated by means of the so called realistic regularization (a suitable mixture of charged 
fermions and bosons compensating the zero-point energy of the vacuum). Е 
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ation starting with the interaction-free waves as the zero- 


order approximations. 
By introducing the equation 


Pen [| K,py (3.227) 
into 
V = y, + ef Куру, (3.22””) 
we obtain an integral equation of the third order 
y= Yo + ef кло as ef Kefky y. (3.23) 


It is welle-known that integral field equations correspond closely to sets of alge- 
braic equations. Thus, (3.23) may possess three independent solutions while only 
one of them is attainable by the usual perturbation treatment (i.e. by iteration). 
It is possible that the presumable extra solutions would account for the great 
variety of mesons not explicable by the traditional tools of field theory. 


Appendix 
The convergence of the sequences y, and p, to a solution y and р may be 
easily proved in the classical case if the integrals (3.17) are all limited (cf. 3.18). 
Let us introduce the quantities 


A, y = Vs — 9 Д, = 9,— Фо | 200.9 

and compute 
| | „у= «к, [4,4 (ш Фо] м. Y+ Yol, (A.2) 
еге E alfa dydzK, ha aM, E 


assuming & >0 for simplicity. Introducing a number N >M we have also 
Me e. ‘age tet ae and Жава, (A4) 
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1 
8 SM (A.8) 
Let us compute now the differences 
Yn „= || x [Pn Pn — Pn—1 Pn-ıl 
к pA (A.9) 


=g [| Кф, — Pa) V, + Pn—1 (Yn — DJ. 


With the aid of (A.4) we have in particular 
f dy K,9o(2) 


| va— yı | Sg?N ТШ! 2К в.) | + | dz 


<enl f dy| [ака + |+ ESE | (Алло) 


< g?N* M + 2g?NM < gN*(gN + 2). 


In the same way 


| vs — vs | SENEN + 2) | [м dz K,(Agy + Yo) | 


+ [а ТЕ К, (Дүр + 90) |) S NN 2) (2: NM + 2M} (Al) 
< g*N*(gN + 2) 28N (gN + 1). 
By repeating the same procedure, we find generally 
| Panti — Pal < ENEN + 2) [2gN(gN + 1)]"—*, (п = 1,2...), (A.12) 


hence a sufficient condition for the convergence of the iteration procedure is 


2eN(gN +1) <1 (A.13) 
ог with N 2М, 
Vas 
pr шы A.14 
& 4M aaa 


КРАТКОЕ СОДЕРЖАНИЕ 


Й. Райский, Регулярная теория поля (квантованного). 

Методом функции Лагранжа, развитым Швинглером, сквантовано нело- 
кальную теорию поля, сконструированную в области заключенной между дву- 
мя гиперповерхностями. 

Обычные канонические законы перестановки применяются к парам 
точек на каждой из этих поверхностей. Не выступают представления связан- 
ные с промежуточными поверхностями. Уравнение Шредингера не существует. 
Обе пограничные гиперповерхности физически отличаются тем, что на них 
происходит содействие с измерительными приборами. Вычислено амплитуды 
вероятности переходов между состояниями, лежащими на этих поверхностях. 
Разобрано некоторый метод непосредственного квантования уравнений поля. 
Подано условия, достаточные для существования аналитических решений 
в случае поля сквантованного, сопряженного с внешним полем. Предложено 
к дискуссии общий случай двух сквантованых полей и ی‎ MU rte He- 
сколько условий необходимых для регулярных решений. 
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RELATIVISTIC QUANTUM DYNAMICS 
OF A SYSTEM OF INTERACTING PARTICLES* 


By J. Rzewuski 


National Mathematical Institute, Wroctaw 
(received July 25, 1953) 


The relativistic problem of interacting particles is quantized by use of an approxim- 
ation method. The solutions of the equations of motion are represented as expansions 
in powers of the interaction constant. In lowest order of approximation they obey the 
equations for non-interacting particles and therefore may be quantized directly. Higher 
order approximations may be expressed in terms of the lowest approximation allowing 
all physically important quantities to be expressed as expansions in powers of the inter- 
action constant with coefficients depending on the unperturbed operators with simple 
commutation rules. Equations determining the S-matrix for scattering problems are 
given. As an example the explicit solution for S is calculated to first approximation. 


Introduction 


Quantum mechanics, so far, was applied either to non-relativistic systems 
of particles, that is systems in which the interaction is of a static type, or to sys- 
tems in which the particles move in a given external field. Extensions to relativ- 
istic systems of particles interacting by means of retarded or advanced forces 
were hindered by the non-local character of the interaction. This non-local char- 
acter manifests itself in the fact that the motion of the particles is described by 
integrodifferential equations and that the energy of the system at a given time ¢ 
depends not only on the values of the dynamical variables at this time but also 
on their values at times ¢’ different from t, the connection of t and t’ depending 
on the special character of interaction assumed in the problem. 

It is an open question whether systems of this type may be put into a canon- 
ical form and then quantized according to the conventional rules. There are hints 
that this is in fact possible (Pauli 1953, Rzewuski 1953), no consequent formalism 
of this type, however, has yet been constructed. | 

In this paper an approximative solution of the problem is proposed, assuming 
the interaction between the particles to be small. The solutions of the integro- 
differential equations are expressed as power series in the (small) coupling con- 
stant. In the lowest order of approximation they are solutions of the problem of 


Nus V oP ee leu Cun Е 
* A short account of this investigation was published in the Nuovo Cimento, 1 
(29) 
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non-interacting particles. This problem is local and may be quantized by the 
conventional methods. Higher order approximations may be calculated, step by 
step, in terms of the lowest approximation. 


1. Variational Principle 


Consider a system of a finite number of particles. Let x, denote the mechanical 
mass of the n-th particle т” Из proper time and q(t”) (и = 1, 2, 3, 4) its world 
line. 

The behaviour of the system is, most conveniently, described bv Fokker’s 
well known variational principle 


OW = 0, (1.1) 
e 1 fs E. 
W = — >" Ү (97) dr" + 5 Уну ед [face 9») 4”, dr dre 
where 
ЖЫ n > 
4, = 2 and (91)? = —1. (1.2) 


In the following units are chosen in which h = c = 1. e, and е, are interaction 
constants. С (q7, q7) describes interaction between the n-th and the m-th particle. 
The summation óver the dummy indices 
is understood. The integration limits 
are denoted by с; and o, which means 
that one has to integrate overt” between 
the values т? and 77 for which the world 
line q(t") intersects two given space- 
like surfaceso, and о». In Fig. 1 the 
vertical axis corresponds to the time 
Фо = 194 = it, the horizontal axis sym- 
bolizes the three space coordinates q}, 


Fig. 1. 7» ds 5 
The four homogeneous coordinates 
4, are not independent, due to relation (1.2). Similarly the four equations of 
motion corresponding to each 4,, which follow from the variational principle 
(1.1) as necessary conditions, are not independent. One of them is a conse- 
quence of the three other. 
For the purpose of quantization it is, therefore, convenient to replace the 
homogeneous variational principle by the corresponding inhomogeneous one. 
This is easily done by introducing into (1.1) the ”three-velocities” 


Р cun 
т uM S, чут e 
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From (1.2) there follow the well-known formulae connecting the four-velocities " 
with the three-velocities v? (greek indices run always from 1 to 4, latin indices 
from L to 3. The convention about summation over dummy indices is adopted 


for both kinds of indices). 


АО ы a e бы 
No ° ¥1— (wp? 
The last two formulae may be written 


n 


v 
d, za) ие (1.5) 
if one remembers that 
nel. (1.6) 
With help of formulae (1.4—6) the variational principle becomes 
OW =0, 


бұ 1 Oz 
W = > «f NIIS (of)? de" t. DA e, en Jet. m di^ di" er 


where again the integration over t” is to be carried out over a region (t, 1), 
where 7 and t are the time coordinates of those points in which the world line 
of the n-th particle q7 = 47 (6) intersects two given space-like surfaces су and оз. 

If the system of particles under consideration is closed (not acted upon by 
external forces) then the variational principle must be invariant with respect to 
translations of space-time. For relativistic reasons it must also be invariant with 
respect to the Lorentz-group. The variational principle (1.1) or (1.7) satisfies both 
conditions if G(g},g”) is a function of the variables 4" and q% by means of the in- 
variant (97— qn. Apart from this condition the form of G is quite arbitrary and 
determines the interaction of the particles; e.g. in the case of electrodynamics 
С (97,97) = à[(qt—q")?]. We shall carry out the calculations for a general С, 
thus including into the considerations also mesodynamics and theories with ex- 
tended particles. 

It is sufficient, without any loss of generality, to treat G(g7,q,) аз a symmetric 
function of its arguments: 


This follows from the possibility of interchanging the integration variables in the 
double integral in (1.1) or (1.7) 


933 ШЕСІ dt^ di" = 


Os (1.9) 
1 | 1 БИ % E 
=F EY een fos toan + ear an 
4 n m Dn 
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It may be noted that the symmetry property (1.8) is not connected with the as- 
sumed invariance with respect to the Lorentz-group. In fact, a function of (9; —4;)' 
may be very well unsymmetric, since the functional dependence on the invariant 
(q? — q”)? may be — without loss of invariance — different within the light 
half-cones of the past and future. 

In the double integral in the variational principle (1.1) or (1.7) we have terms 
with п = m, describing the interaction of the n-th particle with the m-th particle, 
and terms with n = m describing the interaction of one particle with itself. These 
last terms give rise (cf. Wheeler and Feynman 1949 and Feynman 1948) to self- 
energies of the particles which are unobservable quantities and may be considered 
to be contained in the constants x, as mass-renormalization factors. This renor- 
malization of the masses is equivalent to the corresponding procedure in the 
field theory, and allows to write the variational principle (1.7) in the form 


ôW = 0, 


a № „Vi — (vs dir + ; X enen | сату ar de (1.10) 
Та CA n--m сі 


where the constants x, have now the meaning of the experimental masses of the 
particles. We shall use this form (1.10) of the variational principle as the basis 
for our considerations. | 
The variational principle described in this section is of a non-local character 
(even for conventional electrodynamics)’ because of the occurence of double in- 
tegrals connecting particle coordinates (or velocities) at. different times. This 
non-local character, as already mentioned in the introduction, is the source of 
difficulties for a direct quantization by means of the Hamiltonian formalism. 
We shall use, therefore, for the purpose of quantization, approximate methods 
in which the lowest approximation corresponds to a local theory offering no diffi- - 
culties in quantization. | 


2. Equations of Motion 


The necessary conditions for vanishing of the first variation of (1.10) may 
be derived in a similar way as for local theories and have the form of integro- 
differential equations containing time derivatives up to the second order and 
single integrals over the time interval (01,02): 


"A Eee = eu Ff, Eo 


where 


> 


0; 
9А" oA” | 

UE = >. «=_=: АСЕ) . n п жы т m: m m 

dde" sur OD n [ ea md“. — (22) 


т п ка 
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Introducing a new variable T7, we may replace the second order integro- 
differential equations (2.1) by a set of twice as much first order equations 


ж, ve 
т п. 1 
mcr (2.3) 
dt” 
= = 6 FS, (2.4) 


of which the first half are pure differential equations and the second half are 
integro-differential equations. 
Equations (2.4) may be replaced by pure integral equations. Indeed, inte- 
grating (2.4) over the region (¢,”,t”) we obtain 
mi= pl" +e, if wh — aon" Fe de", (2.5) 
бі 
where 
1 for ¢>0 
ret e 
do 0 for t«O. 


Primes denote here dependence on ?"'- рі" is an arbitrary constant representing 


(2.6) 


the value of 7t, at the time #1 corresponding to the hypersurface с. Indeed, in 
view of (2.6), we have 
а Шоты TS (2:7) 
The times t” correspond to a space-like hypersurface о situated in the region 
between c, and o, (Fig. 1). 
We may obtain an equivalent set of integral equations by integration of (2.4) 
over the region (t”, ta”): 


Сз 


mi = pi + e, f qn — m Fede (2.8) 
where: 
0 for £50 
adv un 2: 
7 (2) Ws E 1 99 t < 0 ( 9) 


Неге р? is again an arbitrary constant representing this time the value of лу 
at the time t} corresponding to the hypersurface 6; 
my (t") > pj іс ә 05. (2.10) 
It may be noted that the integral equations (2.5) and (2.8) are of a mixed 
type. They contain integrals with constant limits (cf. 2.2) as in F ош equa- 
tions and integrals with variable limits (447,27) or (t”,t”,) as in Volterra’s equa- 
tions [cf. (2.6) and (2.9)]. | Г 
For further purposes it is convenient to combine equations (2.5) and (2.8). 
We first assume, that (2.5) and (2.8) represent the same solution of (2.4). For 
this it is necessary that the functions лу represented by (2.5) and (2.8) satisfy the 
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same boundary conditions: 
Сз 


PR x 3E 10А. dt” = р. (2.11) 
€, 
We may now take the average of (2.5) and (2.8) getting in this way 
Ф| 
Ере 2l D E(t: =i") ov) Ет а, (2.12) 
0 
where 
1 
Ph = (E + Pe’) (2.13) 
is an arbitrary constant, 
аре" 
= 0, 2.14 
T (2.14) 


and 
Г for T0 
Lora adv == 
e(t) =n" (t) + 7° (t) —l for t<0. 


On the other hand we may subtract (2.8) from (2.5) getting, on account of 
ne (t) — 2d” (t) = Із 


(2.15) 


Ca 


Pi’ — Pe” = е, / Руна КЕ (2.16) 


0 


which is nothing else than the necessary condition (2.11). 


Introducing (2.3) into (2.5) and (2.8), we get 


о» 


in 
n P 5 е, жарға So ret (,п п’ $5 nint Туп! 
Er 1-( + = VIZ | nt m тут; ағ, 
Y (2.17) 
Pk о В 
d deme IUBENT Des REST. 2l "estt yvy Fondue 
n 1 


We may now carry out the second integration 


0% 


1n 

n n р r n n! E Te d 

= qi! + |, (г) Vi — rj de" + 

n 
бі 


TT Оз 
е 2 — ^ ” „ „ " ” 
ite mall ТАН (Ы bat t”) V1 — zn (47 — nc ) s; ER, dt” dt” , 
ao E (2.18) 
qi = ge a = а =) ү1 — (v7)? dt” 


n 
"MET 
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б» 
е, 2 р SS , n ” ” р „ 
Er I N (P у — (qae? t — ew Er de de”. 
n 
бі 


Неге the first equation represents that solution of (2.4) which goes over into q!? 
for t” = t? and the second equation represents the solution going over into q?7 
for t" — t5. Assuming again that both integral equations (2.18) correspond to the 
same solution of (2.4), we get the necessary conditions 


in = S: 
n Pk n n en un "o gp , , , , 
qe = Qf + OE = ffy; =(@ у ет — "yu? FR ааг", 
2: ES (2.19) 
E e ———— 272564) ; 
ar ar [Vi cp ae — [Vy me буурага. 
n e n = 
In order that these conditions are selfconsistent, the following equality must hold 
Pk 3 1 (тұол,т En n\2 ret (n T Vena Tori TIF TT 
z Үт — (vr)2dı + 2 [| vien N (P — "Pv? Fe да 
%, Xx 
о pr (2.20) 


Ge Оо» 
2n iat E E u №. К Ариу 5 
m | Y1— (ar + AE | ү V1 — (0 (^ — P^? Риа" de". 
X Hn 
бі бі 5 


We may easily verify that equation (2.20) is in fact satisfied on account of (2.16) 


and the obvious equation 


те (t) n? (t) —1 (9.21) 
Now, we may express р! and p?” by р" and the quantity 
| ve Ер dt" 
€, 


by means of equations (2.13) and (2.16): 
Oz 
1л On “| n pn dt” 
D Р. v, Py, Ob, 
Tı 


7 (2.22) 
Y e x 
p? = pi^ + a [вуй 
C1 


Introducing (2.22) into (2.19) we get, on account of (2.21) and (2.15), 


O3 
0л = " 
д = а „= f (е-е) У (FJ? а” 
%, 
бі 
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1 , ” ” ” , „ 
е (t Ее ү Ер dt” dt’; 


+ 
X | xD 
—D 
са “ 
а 
~ 
"ts 
28 
— 
I 
[=] 
~ 
м2 
to 
54 


ur (2.23) 
q = q7 + = |! ger — er) Ут гу д” 
^X 


n 
бі 


Оз 
a RECON A 1 , ” ” ” , ” 
б Te CE) V1 — (07)? E E (Б ==) v? Fr dt" dt” 
х 
п A 
The arithmetical mean of ma (2.23) is 


q = д" E — er — t”) V1— (027) dr" 


(2.24) 
Ca 1 1 
+ én = € (t — 6%) Ие (07)? = e(t” — Ра dt", 
EB 2 2 
where 
0л 1 In 2n 
eu) (2.25) 
is a constant of integration. 
Their difference 
ба 
oe — 9h" al V1— (p? ae 
(2.26) 


«уен; e(t" — 1”) v FF, dt” dt” 


is simply the necessary condition (2.19) expressed by means of p?" 
Equations (2.4), in either of the forms we have given to them in this section, 
may be treated by the method of successive approximations. For this purpose 


we assume the solution of these equations in form of expansions in powers of the 
coupling constant: 


ГЭС СЭ (2.27) 


и = 3 4%», (2.28) 
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with 

а" 

ме (2.29) 
Introducing the series (2.27), (2.28) into (2.12) and equating coefficients of 


the same order on both sides of this equation, we get in lowest approximation 


(sn _ 
pon 


TÍO = pi^. (2.30) 
In the first approximation 
(Dn äl n ___ 17) O ря’ тул” 
Us ле; 3 € (t U ) QE : FLU QUA. (2.31) 


а, 
and so оп. Here F{)” is that expression which is obtained from (2.2) by replacing 
qy, q, and 4” s 40» 90” and б The connection between x7 and v; 
in lowest doof pub is from (2.3) 
Я (о)п 
MU 
T” = рр" = (2.32) 


The equations of motion in lowest approximation are from (2.14) 


d sc, 0” 
Ale zo a 


They are just the equations of motion for non-interacting particles (e,— 0) (cf. 2.1). 
Similarly, introducing the expansions (2.27,28) into (2.24), we obtain in 
]owest approximation 


On f^ 
d an PE | кб n ER an 
x 


n 
б, 


T 
qm + v9» [e EY (tt + 4)] _ (2.34) 


and in the first approximation 


б, 


б )س‎ 
dior |е VIET ars 
x 


n 


(2.35) 


O1 


EE se (feet) Vi Come — g (p! — m") on" FOR” qe de”, 


where ۷1 — (47% — (0 The Md loud the first order term in an expansion of the root in 
powers of the interaction constant. 
The equations of motion in lowest approximation are, on account of 
j 


> e (t) = 0(0), (2.36) 
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КЛ ee 


dq (0 

dt” 

which is equivalent to (3.32). 
Instead of expanding equations (2.12) and (2.24), we can as well expand 
equations (2.5), (2.8) and (2.18) getting in the lowest approximation equations 
of the same form as (2.32), (2.33) and (2.37). The upper index 0, 1 or 2 denoting 
constants of integration (which in case of 1 and 2 have the meaning of boundary 
values of the corresponding quantity at o, or сз resp.) must not be confused with 
the upper index (0), (1) or (2) in (2.27—37), which denotes the order of approx- 


imation. 


ры TEQUE 1— (92, (2.37) 


3. Quantization 


The equations of motion in lowest approximation are, for all types of solu- 
tions described in the preceding section, of the same form (2.32), (2.33). They 
are derivable from the variational principle 


SFO - s jv — (y? dt” (3.1) 


as necessary conditions for the vanishing of the first variation. This variational 
principle has local character and allows, therefore, canonical formulation and 
quantization by the usual methods of quantum dynamics. Since there is no inter- 
.action, the motion of each particle is independent of the motions of other particles. 
Тһе system may be described by a number of canonically conjugate variables. 
The canonical coordinates are q? (i = 1, 2, 3), the corresponding canonical mo- 
menta are 


p; EUREN a: f 3.2 

۲ (3.2) 
The equations of motion (cf. 2.33) have the form 
dp? и 

cmn 3.3 

= (3.3) 


(for simplicity we leave out, for the moment, the upper index (0)). 
The commutation relations between the canonical variables may be derived 
by any of the conventional methods. They are 


[92,98] = [p?, pe] = 0, [FPF] = 50" д. (3.4) 
These commutation rules hold; at first, only if the arguments t” and 477 correspond 
to the same space-like hypersurface о. The commutators for arbitrary times follow 
from (3.4) and the equations of motion (3.2), (3. 3). 
From the constancy of p? we deduce ШЫЛ мул that the commutation. 
rules containing pz hold unchanged for arbitrary times: 


2, ра) 50, [47, ри = io б, (for arbitrary times) (3.5) 
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To calculate the commutator (47, q7] for arbitrary times we have to express 


one of the quantities (say 9”) by its initial values at a time corresponding to the 
same с as the time t". From (3.2) and (3.3) we have first 


n 


dq; 
dt” 


a 2 Тү кене 
рь = const., v, = 


= const, gi (t") = а (0) + oF (Г — P). (3.6) 


In the solution (3.6) we may express v? by p? by means of the equation (3.2) 


DAS pou ; 
Vas + (PI)? 
It may be noted, that the inversion (3.7) of (3.2) is true also for operators, since 
the v; commute with each other and with the p? on account of (3.5). 
Introducing now (3.6) into the commutator [4”, q7] for arbitrary times and 
using the relations (3.4) and (3.5), we may write 


[z (P). FE] = [aR (9), ag (7) + pe —m = 
Ouen] с? (3.8) 


To evaluate the commutator in (3.8) we may use the relation (3.7) and the well- 
known equation | 


(3.7) 


д т 
4. Ар = 1 шеш (3.9) 
Pi 


following from (3.5). We get so the result 


j p | Шү 
| رھ‎ a |= — V1 — @)? (6, — 00) (8.10) 
ЕЕ ale a+ (P| x TAS 


and 


— m Pi Pk пут 
с sperme e 
[92 ( ), ак (t”’)] al E эс dl 


(3.11) 
= ү1--( (a — 99 (7 — 2"). 
x ` 


Having thus established the commutation rules in lowest approximātion (no 
interaction), we may proceed to the calculation of the physically important tran- 
sition amplitudes between states on 01 and 0з. For this end we apply here — for 
quantum dynamics — a method developed by Yang and Feldman (1950) for the 
purposes of the field theory. j 


The canonical commutation laws (3.4) are satisfied by each of the independent 


2n 2 : E 
sets of canonical variables: 40, р”, 4", рі", and Фар. According to a well- . 
known theorem (Dirac 1947) they must be connected by a unitary transforma- 


tion. We may write, therefore, 


z = 3 à 
pP = Spy” 5, sts =1. (3.12) 
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The unitary operator S transforming operators on c, into operators on 02 
may be interpreted as Heisenberg’s S-matrix. It transforms the corresponding 
state vectors (с) and (сз) of the two unperturbed systems according to the 
equation 

Фо) = 514 (01) 0 (3.13) 

and Из matrix elements give the transitions amplitudes for transitions which 
occur between these states. 

To calculate S explicitly we combine the first equation (3.12) with condition 
(2.16) 


95 
ірі: 6 Se, | v: Fp ар, 
Cı 


у, (3.14) 
[pr е ДЕЛЕ 
Taking the arithmetic mean of these two equations we get 
g? 
1 : 
LEES E ri | v FS an (3.15) 
where 
[4, В| = AB —BA, {A,B} = AB + BA. (3.16) 


A second equation for S is obtained analogously from the second equation 
(3.12) and the necessary condition (2.26): 


lax", 5] = 
1 n 9% Оз i 
=; Ís, et: if Vi — (аа || |: У — =" и) v Fi a ae" 
Ha Xn 


(3.17) 
Equations (3.15) and (3.17) determine the matrix S uniquely. 

The considerations of this and the foregoing sections enable us to express 
all physically important quantities as expansions in powers of e, with coefficients 
depending only on the local quantities 429”, p(9", 009" satisfying the commutation 
rules (3.4), (3.5), (3.10) and (3.11). Thus the relativistic problem of interacting 
particles seems to be solved in principle in the frame work of quantum dynamics, 
without use of the field concept. It may be regarded as an extension to quantum 


theory of the ideas of action at a distance advocated by Wheeler and Feynman 
(1945, 1949). 


4. Applications 


To illustrate the general considerations of the foregoing sections we shall 
calculate the S-matrix up to the first approximation. For this purpose we expand S 
in powers of the interaction constants e,: 
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5 = № SS (4.1) 
s=0 


and introduce (4.1) and expansions (2.27, 28) for 4; and v? into (3.15) and (3.17). 
Equating coefficients of equal order in the interaction constant on both sides of 
these equations, we get in the lowest order 


[ру оо =; (4.2) 


о» 


ы 1 p% pi SEP Tri img 1 
fa, 5) = 2 |89, a | Var) - 2159. | {таг (аз) 


бі 
and іп the first order 


1 f 
ірі”, SO] = 9 ів, е, ТА ӘЛЕН rar) (4.4) 


бі 


1 On P: 1 On % 
(а до] = a ізо, e [vi — (vr)? am db зе | ү: — (we ar 
n Xn 


Oy 
1 е — 1 e у ы ; 
un (5%, an | Vi — Gf er г) FO" д). (4.5) 
Hn 


From equation (4.2) we infer that S® can not depend on 41%; 


Эр" | 


SO = SO. . C ERA 
Introducing (4.6) into (4.3) one obtains a set of equations determining the func 
tional dependence of S® on р?” (cf. 3.9) 
: Q SO 0m Or P 
Гай”, 5% [p] Peace PT E Sr | VI de. (4.7) 
х 


` The order of the factors on the right is arbitrary on account of the permutability 
of ро” and v”. 

Expressing now the right hand side of (4.7) in terms of p?" with help of the 
equation 


Vs (тра ы ы едын (4.8) 
Vz + (py? 


following from (3.2,7), one gets | 
: on Os о» ER, 2 Ca 
Pk ü Pk kr. ү? Ony2 jn 
+ | y1 — (v2 de? = | а = a x» + (p;”)? dt”. (4.9) 
ale er ит: 
с бі = €, 
Equations (4.7) take now the form 


Os я 
(0) | 
4 95 44 secs | ж E (p)? dt” (4.10) 
др?! Op; 
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Integration yields 
ОЕ nA И Fe 
go. CZ JV Gran (4.11) 


m Cı 


The additive and multiplicative constants of integration are chosen in such a way 
as to give 
[S =, (4.12) 
corresponding to the statistical interpretation of the S-matrix in lowest approxi- 
mation. The initial state remains unchanged in the course of time due to the lack 
of interaction. 
To obtain the first approximation we introduce (4.11) into equation (4.4) 


о» 


pte 2 ро? 2 dtm 
Ірі», 59] = iz j VA оа жа 1 000" FOR л . (4.13) 
2 


т о, 
о, 

To simplify the problem, we may now make the assumption that the surfaces 
сі and c, are situated in the remote past or future respectively and that the col- 
liding particles are in large distances from each other at times corresponding to o, 
and сз. In this case the interaction of the particles at су and с, is negligible, which 
may be used to simplify the classical equation (2.16). Indeed, the right hand side 
of this equation is, according to (2.2), 


C» 


ЖЕНЕ 
T aie [Ar E= Я qe 
$ да? i 


mn 04 


Here the second term contains a total derivative with respect to time and may 
n 
t2 t2 


be integrated to give 
„| G(T — Ф") 0") = —e, A (GD, . . (415 
zz > 2 б ü | ti 
According to our assumption about the interaction of particles at сү and А this 
term is negligible and (4.14) becomes 


n n гы 9G rg) EM n m 
JE Fr di =e, Уе SIE ааа. (416) 


man 'g 
Applying now the quantization procedure-to the simplified equation (2.16), we 
get in place of (4.13) the equation 


(4.14) 


о; 


On 
£ 1 --іу [\х? + (p%™)2 ат д 
[222,59 = 5 (е 2] MORALE) > aq ^ КОО ‚ (4.17) 
5 k д 


04 
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which is satisfied by 


. is s 2 Om 2 дрп 1 P 

uet i EVE үи ) n 7 

2 | та 32,4 | Аат + утро] | (4.18) 
n p ў 


The additive function of the momenta ЛР] occuring in (4.18) is determined 
by equation (4.5). 
Collecting now (4.11) and (4.18) and making use of (2.2), we get finally the 


following formula for the S matrix up to the first approximation 


sa _ 216 -іх [VAL ge po)? des 
D 


m o, 2 


(4.19) 


УУ Enem Í ТЕ R арс отш T +. 


mom 


This formula. differs from the corresponding formula (16) of the note published 
by the author elsewhere (1953) as the exponential factor of modulus unity is there 


missing. 


КРАТКОЕ СОДЕРЖАНИЕ 


Я. Жевуски, Релятивистическая, квантовая динамика системы вааи 
модействующих частиц. 

Методом приближений сквантовано за дачу. релятивистически взаимо- 
действующих частиц. Решение уравнений движения представлено как раз- 
витие постоянной взаимодействия на ряд степеней. В первом приближении 
исполняют они уравнения для не взаимодействующих частиц и поэтому 
могут быть сквантованы непосредственно. Высшие приближения можно 
выразить при помощи выражений найнизшей степени; позволит это все, 
физически важные величины представить как развитие на ряды степеней 
постоянной взаимодействия с коэфициентами зависимыми от не возмущен- 
ных. операторов C простыми правилами перестановки. Представлено ма- 
трицы $ для проблемов рассеивания. Для примера вычислено непосред- 
ственно решение для 5 с точностью до первого приближения: 
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NOTE ON THE PROBLEM OF COORDINATE CONDITIONS AND 
EQUATIONS OF MOTION IN GENERAL RELATIVITY THEORY 


By Roman TEISSEYRE 


Institute of Theoretical Physics, University of Warsaw, Warszawa 


(received September 2, 1953) 


It is shown that in the approximation method of Einstein and Infeld (Canad. J. Math., 

1, 209, 1949) the equations of motion in newtonian approximation and in the next one 

are independent of an arbitrary function a, appearing in the solutions of the field equa- 
3 


tions. This result together with the results of Einstein and Infeld concerning a, and bọ 
4 5 
lead to the conclusion that in their new approximation method the equations of motion 


in newtonian and first beyond newtonian approximation are uniquely determined by 
adopting the newtonian potential as solution of the first field equation. 


In 1949, Einstein and Infeld have presented in a new, rigorous form their 
method of deriving the equations of motion of particles from the equations of the 
gravitational field (Einstein and Infeld 1949). It was this work which was the 
starting point of the present investigation. The same notations as there are used 
here. 

It is well-known that in every step of the new approximation method for the 
field equations there appear some arbitrary functions; their appearance is connected 
with the existence of the Bianchi identities. Einstein and Infeld (1949) discuss 
the meaning of these functions and their bearing on the equations of motion. 
More thoroughly this problem is treated in a paper by Infeld and Scheidegger 
(1951) and recently by Infeld (1953). These arbitrary functions correspond to 
some arbitrariness in the choice of the coordinate system. 

In the approximation method of Einstein and Infeld the solution of the field 


equations in the (4-4) th step will be given by the functions Yoo Уот» Ут" Етот 


these solution we can build a more general solution in which the equations will 


be satisfied by the functions 


Ж. T Yoo% (la) 
k—2 k— 2 
yd i Yom zu ғо т, (1Ь) 
k—1 k—1 k—1 
/ = я — д a.. + д а» (1с) 
} mn Ўт = imn i inm mn Rot ih 


where a,, a, are arbitary functions. 
kel К 
(45) 
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Similarly for the set of solutions y,,, Усть Ymn We can construct more general 


k k+1 k 
solutions, by taking 
ТА а Yoo + b, o (2a) 
k k k 
: y = рот dus UTE (2b) 
k+1 k+1 k+1 Ю 1 
ios = Утп zs boss zi OFF, FIT Sinn b, s (2c) 
k k k k k 


We note that we can always put а, = b,, because the solution y',,,, = Ymn + amn t 
k k k k k 
“ға, m— Omn а, changes the field equations and induces new solutions for у, 
k’ ke k 


Yom Which are analogous to (2a) and (2b). 
k+1 


Einstein and Infeld have shown that the equations of motion in newtonian 


approximation and in the next one are independent of a,, b,. Their proof rests 
| ; 4 5 
on the assumption of independence of the integrals 


$ A uk ni day t= 475, 6 


of a, b, They have also extended their proof to the independence of a, b 


4 5 | 21—2 21-1 
up to the (/-1)th order of approximation. 
In connection with Infeld's work (1953) the question arises as to the anal- 


ogous independence of a,. This independence is proved in the present paper 
5 : 
by direct computation. The first step consists in computing 24,,, and the field 
6 
equations to the fifth order of approximation without recourse to coordinate 


conditions. The introduction of these conditions would simplify the calculation 
but, on the other hand, it woud be equivalent with assuming some definite func- 
tions for a,, a,, b,. The second step gives the proof of the independence. 


I 


Our calculation of 2 ^,,,, leads to the same expressions as in the paper of Ein- 
4 ; 


stein and Infeld (1949, Eq. (A. 12,3)) except for the addition of 


(ауы ml as) = ( “т Уһ, a Ф nY ms,s —29 Vas) d 
% 5 1 4 2 4 4 4 


x (2, 20 Yom;n E 29 o Yon,m Fw Yos,s Уот,л > Yos,s Von) y 
k 13 13 -3 3 3 3 


T (yt Voss — 2 eor 
3 3 13 
where y,, = 29 
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All these additional terms vanish by the Einstein-Infeld coordinate conditions 


Yms,s = 0, Уов eS Yoo о” 
x kic kid 
Analogously the additional terms for 24%, (i = 4,5) are 
2 lor : 0 
4 
24. 0 
4 
2 en OMA Yos,s vs 2g oP m 
3 3 
П 
In order to show the independence of the equations of motion of a, it is only 
necessary to prove this independence for the surface integrals З 
$ A, n, dao. (i = 4, 5, б) (4) 


To prove this, we put into the integrand instead of у, Yan Ymo the more general 
JOUÉ DIETE 


solutions 


D 2 , oy @ бы 
О ув у om Yom s аот?» У mn = Ymn n о а оо" (5) 
2 2 3 3 3 4 4 3 Ж 


This changes also the differential equations for y,,, Yomi for Yon, we get for in- 
4 5 4 


stance 


, — 
Y; 00,55 = Уоо,55 di Qo oss 
4 e) s! 


and from the arbitrariness of a, it follows that 
3 


RS = Yoo ne ao o (6) 
4 4 27,1 
together with (5) will form the solution of the field équations. The knowledge 


of Yom is not necessary. Starting from the differential equation 
5 ч 
Vo sm Yom,ss mus (Ф а) Pu (Pm UR (7) 
5 5 3 3 
we get in the same way for Ут. 
5 


D 
y om, 77 Yom,s "s Ps som Фат ао, ei Eskr Cnr k* (8) 
5 5 


Here the last term arises from integration of the divergence (7) but it gives no 


contribution to our surface integrals. 
As shown in the appendix, that part of our integrands which depend on the 


arbitrary function a, éan be put into the form F4, , where Ри, = Ё, Then, 


3 * 
it is easy to show (Einstein and Infeld 1949) that 
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$ Бат, D day 7 = 0 (9) 
From this the independence of a, follows. This result together with th above 
3 
mentioned result concerning a, , b, can be expressed as follows. “ 
4 5 


The equations of motion up to the 2/ th order of approximation, i.e. the 
integrals 


$ uk Пк Kay? for i < 2l, (10) 


do not depend on a,, a,, b, and therefore do not depend on the coordinate trans- 
21-3 21-2 21-1 
formation involved. 


This generalization rests on the fact that a,, a,, b, enter the integrals (10) 
2 3 2-2 21-1 
in the same way as a,, a,, b, enter the integrals 
СИ М 


$ nd,‘ for i <6. 
i 


Thus we see that in the new approximation method the equations of motion 
in newtonian and first beyond newtonian approximation are uniquely determined 
by adopting the newtonian potential as solution of the first field equation. 


Appendix 


The part of the integrands which depend on the arbitrary function can be 


put into the form F,,,, where F,,, =—F,,. For 2 „р we obtain 
6 


(49 os a, бк сақ Ag ok 420,2), зщ (2ф, Qo Omak қ 20, ао “ЖАР 
13 13 22d 3 1 
PE (2, йок д > 2%, ао; Omk),s T (Yok a, sm — Yos а), 
13 13 3 3 SIS z 
ат (Yor ao rk Ona — Yor Чо, к дик) з 55 (vo, Qo rr Omak — Уо Qo rr РЯ 
3 3 3 3 3 3 3 3 
J (Yost Ч от — Yok,s а дк = ее Cok — Yomk ME 
3 3 3 3 3 3 3 3 
+ Е йок даз — Yos,r o,r 0 357 4; Ж а; Sak — Yor,r %o,k DE 
3 3 3 3 3 3 3 3 
F (2765 mk Ar 2Y ok ms%o),s = (aox ao, sm — Qo, а im) s 
3 3 3 3 3 3 3 3 
y" (а, Q, rk Ong IX а, а. rsOmk),s as (a, rr ао бок RE Qo rr Сок б) : 
3 3 3 3 3 3 3 3 
Similarly for24,,, we get 
4 


(a, ok ms > a3 os mk) ‚5? 
81 31 


and for 24, 
5 


(o, Cok Fel Фк а, )) „» 
3 3 
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КРАТКОЕ СОДЕРЖАНИЕ 


Р. Тессер, Замечание к вопросу условий для координат и уравнений 
движения в общей теории относительности. 


Доказано, что в методе приближений Эйнштейна и Инфельда (Canad. J. 
Math., 1, 209 (1949)) уравнения движения в приближении Ньхтона и вследую- 


щем приближении — независимы от произвольной функции ар, появляющейся 
3 
B решениях уравнений поля. Этот результат совместно C результатами Эйн- 
штейна и Инфельда относящимися к Ag W by приводит к выводу, что в их новом 
4 5 
методе приближений, уравнения движения в приближении Ньютона и сле- 
дующем, определены однозначно признанием потенциала Ньютона как pe- 


шения первого уравнения поля. 
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An equation for the eigenvalues of the energy of two body systems is derived from 
field theoretical first principles. The simple model of non-relativistic spinless charged 
fields is treated as an example. Since the energy eigenvalue problem refers to a fixed 
time instant, the explicit consideration of retardation may be avoided and the instan- 
taneous electromagnetic field may be represented as a sum of two terms: the first in the 
form of free waves, the second in the form of a stationary solution. The energy field 
operator splits into several parts: the kinetic energy of charged fields, the energy of 
photons, the interaction energy between photons and the charged fields, several self- 
energy terms (due to the interaction between a charged field and the electromagnetic 
field produced by the same charged field), cross terms (due to the interaction between 
a charged field and the stationary electromagnetic field produced by the other charged 
field). The cross terms may be transcribed to the configuration space and yield, besides 
the Coulomb interaction and the Breit terms, some other corrections. The equation 
for the energy eigenvalues, i.e. the generalized time independent Schródinger equation, 
is written in terms of relative coordinates of the two particles. 


1. A simple field model 


In this preliminary paper we are dealing, for simplicity, with non-rela- 
tivistic, spinless charged fields үу, v4 coupled to the electromagnetic field A,. 
We may assume the following! Lagrangian 


2 
1 cise 
"e » Т es gs. rad — e, Avi (багай + e, A)v, 


д 1 E д ж 
a ET ( EA. е) | 3h m «oid bs L, т. (1.1) 


2 dt 
with Кр 
дА 
L, er A^ = ч (EL) 


Ox, Ox, 


The field quantities y, and A, are functions of X and t. We may consider 
V1 (фә) as the electron (proton) field. 4, is the electromegnetic field (41, Ag, As) = 
Mac MI ا‎ «у ا ر‎ LLL ee на 
1 Einsteins summation convention, and units ћ = c = 1, are used. 


(51) 
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The generalized momenta for the charged fields are 


IL De wow IL i 3 
S reese ESS , Т, = VD aeri n AR: a? where Ls d el: 122 
T 30823 Par Ta 2 f (1.2) 


a 


The charge and current densities are given by 


2 2 
pres i M e, (у T,— Ta Ya) FF 3 a V. Vs (1.3) 
a=1 а=1 
2 
б = V C [чу (i grad + e, A) ya + (i grad —е,А) v; va): (1.3) 
Ј Ж Әт, D ( gra a Wa a a a 


In spite of the fact that in this non relativistic field model the charge and 
current densities do not form a four-vector, we shall use occasionaly the notation J, 
with и = 1, ... 4, where /, = ie. The charge and current densities (3) satisfy the 
continuity equation 


9,], = div] +e — 0. (1.4) 
The energy of the system of fields is given by the Hamiltonian 
Hs f @xH, (1.5) 


where 


1 br 
В =) = әт (grad —e, A) v, (i grad + e, A) v + e, ey y Hem. (157 


7 
with 


mes | 
Has = (4,4, + grad А, grad A,). (1.5”) 


The momentum of the charged fields is 


P= -X fern. grad y, + grad у, т.) = 
2 (1.6) 


i * 
ы ақ endet 34. ен. ТАС 
E ) fe x(p, grad Ya — grad y, > Ya). 


The field equations derivable from the Lagrangian (1) by varying the field quan- 
tities y, are the well known non-relativistic Schrédinger equations 


1 
әл {grad + еА) v, + egy, = iy, а-1,2, (1.7) 


while for the electromagnetic potentials the field equations are 
L4, =—/» и = 1...4. | (1.8) 


The above described fields may be subjected to field quantization. We assume the 
usual commutation relations for the electromagnetic field and the following 
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commutation relations” for the charged fields 


[va (ж, 0, ү» (X°, 9] = д. (x х), [pa (00), 4х9] 
= [va (е0), v; (x, )] = 0, (1.9) 
whereby the field equations (7) agree with the equations 
ij, = [vo H]. (1.10) 


It may be easily seen that the same field equations (10) or (7) are also obtained 
if we replace in (9) the commutators by anticommutators. In a non-relativistic 
field theory both methods of quantization corresponding to the Bose-Einstein 
or to the Fermi-Dirac statistics are possible. 


2. The Ё and x representations 


The momentum representation (%-representation) is most extensively used 
in quantum field theoretical considerations. To introduce the R-representation one 
used to perform a Fourier transformation 


р (x, t) = олу ? dk y (k, t) e'^* (2.1) 


It is immediately seen that the commutation (or anticommutation) relations for 
the Fourier transforms are exactly the same as those for the field quantities them- 
selves: 


[va (k, t), v; (k’, 0], = д, 99 (x — x), [v (х, 0), рь c, 0]. 
= [v; (х0), v; (x, 2], = 0. (2.2) 


The operator of the total charge is 


Er Ds Ea f d?a y, (X, ї) Ya (¥; 0 = p е, f Pk y, (k, t) palk, t), (0.3) 


while the operator of the momentum (1.6) is 
P = № f dk k у’ (k, t) v, (k, t). (2.4) 


The physical interpretation is facilitated by the introduction of a ” periodicity 
box” in the x-space. Then the functions y (k, t) change into sets Y(t) satisfying 


2 The fact that (9) guarantee the consistency of (10) and (7) justifies the choice of the commu- 
tation relations id It is interesting to notice that in our case we have, from (9) and (2), [ya (х, 2). 


Е ^ ich di — l "canonical" commuta- 
ль (x^, 0] = E Ó;50 (x — x’), which differs by a factor 5 un the usual "ca 
: і if d a non-symmetrical 
tion relations. It is immediately seen that the factor auia disappear if we use y 
# 


Lagrangian with a term yj" (: SET “4 ”) y replacing the two similar terms in (1). 
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the commutation relations 
(Уы (4), V y(t) J. = даь бу. (2.5) 


The operators of charge and momentum become 


Q = Уз 2 % м (д, P= X k, Nf? (t), (2.6) 


М (4) = V, hi (t) * Pa ri (t)- (2.7) 
As is well-known (2.5) may be realized (at a fixed though arbitrary time instant tọ) 
by the well-known annihilation and creation matrices, whereby №“) become 


where 


diagonal matrices with integers 0, 1,2... in the case of commutators, and 0,1 in 
the case of anticommutators on the diagonal. This makes possible a corpuscular 
interpretation of the quantum field formalism. The eigenvalues of N (® (to) mean 
the numbers of particles associated with the field with the charge constant e, , 
mass constant m,, and momentum A; at the time instant tọ. Since we are working 
here in the Heisenberg picture, the operators N(® are time dependent. The states 
of the field may be described by means of the eigenvalues N(?' ofthe operators N(®: 


| NOR, МӘ”, es > (2.8) 


The above ket vector denotes a state of the fields for which a measurement (per- 
formed at tọ) yields the information that N’ particles possess the momentum k4, 
КӨ” particles possess the momentum Еҙ, etc. At the same time the total number 


of particles is known to be > (N@’ + МӨ” + ...). As we are practically interest- 


ed only in states for which nearly all degrees of freedom (Е,) are unoccupied, 
it is convenient to write down in (2.8) only the eigenvalues different from zero. Thus, 
LED and rk) > (2.9) 
denote states in which only the n-th degree of freedom is occupied by a single 
particle either of the type a = 1, or a = 2. (Particles of different types are sepa- 
rated by a hyphen). In the same way 
SRS ROR si, ARREST ЕЗ» (2.10) 
denote two particle states in which either both particles of the type a = 1 possess 
equal momenta £, , or both of the type a = 1 possess different momenta Ё„ and R, 
or a particle of the type a = 1 has the momentum 4, while a particle of the type 
а = 2 possesses the momentum &,,. An arbitrary state vector may be represented 
as a linear combination of the normalized vectors of the above types 
| >= ak; | Aas > ran; |26; Do ee 
+а | Ras > + of: CES 
ырк | (2.11) 
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INKS INE SAN. ^ ~ 
+ ар, Kj os; Re A NU EN bU АМЕ ect a 


The representative of the vector |) 


«МЕ, МЕ... | У = al INS (2.12) 


ик 

is the probability amplitude for finding in a measurement performed at tọ the 
corresponding numbers of particles. In (2.12) and in the following we assume 
only one type of particles for simplicity. If | > is a linear combination of one par- 
ticle states, then only the amplitudes а are different from zero. If we remove 
the box of periodicity, then the set a}; becomes a function a(k). In the same way 
if | > is a linear combination of two particle states then the amplitudes which 
are different from zero are а", or a(R, l). Since the order of the indices A; is 
immaterial, the function a!!(k, J) must be symmetrical (or antisymmetrical in the 
case of a quantization with the exclusion principle). 

In exactly the same way the x-representation may be introduced. While the 
R-representation is connected with the eigenvalues of the momentum, the x-re- 
presentation is connected with the eigenvalues of the operator 


х@= У, Г Фик, У f а.” (k, 1) grad v, (k, 1) 


— grad y,” (k, t) · v, (k, t)]. (2.13) 
Denoting l 
NO (x, t) = pa” (X, t) pa (x0). (2.14) 


we may write 
Q z 3 ЯШ (x,t); X (t) -Y jo NO (x, t) (2.15) 


which is similar to (2.6). The physical interpretation would be facilitated also in 
this case by introducing a periodicity box" (in the R-space), and we could re- 
peat the above considerations by replacing everywhere Ё by x denoting the (dis- 
crete) positions which may be assumed by the particles. In the one particle states 
only the coordinates al are different from zero. When removing the auxiliary 
box, the set ау, becomes a function a!(x) which is the well known Schrödinger 
function for one particle states. The two particle states are described by 255 
or all(x, у) which must be either symmetric or antisymmetric in the case of iden- 
tical partiCles. It is easily seen that the matrix of the total charge (the total number 
of particles of each type (is diagonal together with the matrix of position X given 
by (2.15) or together with the matrix of the total momentum P. This is true because 
the eigenstates of X and the eigenstates of P are at the same time eigenstates of 
the numbers of particles of each type а = 1, 2. P and X both commute with Q 


but do not commute with each other: 


[X; (0), P) = id, / d*xy* (x, t)y(x, t) = ij 3 М T 3 210) 


-- 
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The states of a given total number of particles of each Тере ranir inte Sates 
of the same number of particles when transfering белі the = t the Amp 
sentation. In particular, the ome particle veeter | le >) а бонат eomhizatien ef 
one particle vectors im the &-repres i | 


а> = fd па>дь д em 
Let us compute the transformatiom matrix 
ғ. x) = Зар em 


We have 
РП а> =| b> аја рә Па ey 
om account of the fact that the eigenvectors of more tham eme particle states (Gm 
either of the spaces x or А) are all exthegenal te | 1 Ах Qn. the ether amd ber 
integrating (1.0) by parts we get 
Р а) eye, дагъ у^, ау 
= — dc y^, 9 grad" pe’, 9 | e (x | 1a) 
= if Bx’ f dwp* (x^, 2) grad” Fk lL ақ 
АСТУ rî |= | 
| = —i е |1) grad” Съ | I 
where the properties of y (or y*) as annihilation (ar creation) operators have heen 
wed. By multiplying (218) and (219) by Che] we find = 
ачаа аса e. 


Energy of Bound States 57 


which may be generalized for an arbitrary function of x: 
(Фар а, ) fe), t) | 1x» = fe) < | lx. (2.24) 


Similar computations yield 


(| J x’ grad’ уж” fix yla’, t) | lx = f(x) grad < | 1x > (2.25) 
2 | f d3x'v*(x*, ӘК? grad y(x’, t) | lx > = — grad f(x) < | 1x X. 


3. The transition to the configuration space 
in the case of an external electromagnetic field 


Let us assume in this section the electromagnetic field is an external field. 
In this case the electromagnetic part L,,, in (1.1) and H,,, in (1.5) disappear. 
The equation for the eigenvalues of the energy becomes 


Н |> = WI (3.1) 
fot the ket vector or 


E т (3.1”) 


for the bra vector. Let us assume that at a given time instant the system consists 
of one particle of each type a = 1, 2. In this case only the amplitudes < | 1x,; 
lx,» or <1x,; lx, | > are different from zero. By multiplying (3.1’) to the right 
by | 1x,; 1х,> we get 


1 А " 
fe} | = 3 ered —e, ду; (i grad + е, А) ү, 


t (3.2) 
te А | 1x1; 1x55 = И < | Пхи; 155 
and, by applying formulae (2.22), (2.24), and (2.25) (whose generalization for two 
distinct fields is immediate) we get 


= ( grad, — e4 (ху, 8)? + сүр ER (3) 


+ Е (i grad, — eg A (xs, t))? + canes) | «| lr; lx» = W< |1x,; 1x1) 


2m, 

which is the well known Schródinger equation for two charged non-relativistic 
particles. It is immediately seen that the representative < | 1x,; 1x5? is the complex 

conjugate of the Schrédinger function 
< | Ley; 155 = v* (xi; хз) (3.4) 
while " 
(1x; lx, | >= ф(х, хз). (3.4) 
In this way the connection between quantum field theory and quantum 
mechanics in. configuration space is established. It should be stressed tha tthe rep- 
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resentation applies to an arbitrary but fixed time instant tọ. If the electromagnetic 
potentials are time independent then we may forget also the time dependence of 
the representative < | 1x,; lx,>. But if the electromagnetic potentials are time 
dependent then the state of the system will also change with time and the represen- 
tatives of the state vector may change. It may happen that in the course of time the 
probability amplitudes for finding more than one particle of each type will become 
different from zero. This means the possibility of creation of pairs of charged par- 
ticles by a time dependent external field. Thus, the energy eigenvalue problem is 
meaningful only if we consider the electromagnetic potentials at the same time 
instant & to which the representative < | 1x4; 1x35 refers. 


4. The eigenvalue problem in the case 
of a general electromagnetic field 


The problem of the eigenvalues of the energy in the general case of two charged : 
fields (e.g. the proton and the electron fields) coupled to each other through the 
interaction with the general electromagnetic field has been investigated years 
ago by Breit (1929) and later by several authors who investigated further correc- 
tions to the formulae of Breit. Howeyer, none of the authors presents the eigen- 
value problem as a whole from a unified point of view. It seems to be worth while 
to derive the eigen-equation for the energy from field theoretical first principles, 
to interprete the various terms, and to solve the eigen-equations to a desired 
degree of approximation. This is the aim of the intended series of papers. This 
first publication deals with the simple case of spin-less non-relativistic charged 
particles coupled with each other through the general electromagnetic field. The 
essential features of this new presentation (in particular a trick making it possible 
to avoid an explicit use of retardation) is applicable also in the case of relativistic 
particles obeying the Dirac equation. 

The energy density (1.5’) may be writtten in the form 


2 
HH uo eee Ad а №: = Pa Pa (4.1) 
or > 
H = Hy, + H, m —A,j EIAS (4.17) 
; xt „+5 (4. 
where 
Hy = ух - grad y; grad y (4.2) 
4 2т, 5 = k 
and 1 
dime 
k= odds TM (4.3) 


The electromagnetic potentials satisfy equations (1.8) which may be solved in 
several ways. For our purposes it will be particularly convenient to write the 80- 
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lution? of eq. (1.8) as follows 
А, (x) = A(x) + AÇ (x) dr TD) Aa), (4.4) 


where x denotes in short the four variables x, t. A® is a general solution of the 
wave equation 


04% — 0, (4.5) 


A is a special solution of the Poisson equation 


А AC? (x, t) = —j, (x, t), (4.6) 
where £ is to be regarded as a parameter. The function D^(x, x’) is the Green func- 
tion given by 

== er TS E 
Dx) ED (x — x) for i» tovt (4.7) 


0 elsewhere, 


where D(x) is the well known Jordan-Pauli function. D^(x, x’) satisfies the equation 


OD (x, x") = — P(x — x’) (4.8) 
and vanishes together with the time derivative for t = tọ: 
, 9 + 
D(x, x"), iq = Р(х) = 0. (4.9) 


Indeed, with the aid of (4.5) —(4.9) it is easily verified that (4.4) represents а gen- 
eral solution of equation (1.8). Moreover, up to the first time derivative the so- 
lution reduces for 2 = ty to 


А, (x; te) = AD (x, to) + Ae (х, to), 
A, (x, %) = Ara to) + ИЕ A 
The special solution of the Poisson equation may be assumed in the form 


1, (к, 0 
| — x| 


(4.10) 


A(x, t) = = аза (4.11) 

Solution (4.4) involves retarded potentials for t > tọ and advanced potentials 
for tg >t but if we are interested in the particular time instant / (which is just 
the case for the energy eigenvalue problem) the effects of retardation do not 
matter and we may use the solution in the form (4.10) without any loss of gen- 
erality. I tseems that this simple fact was not seen clearly enough by previous au- 
thors. At first sight one could object against this formulation because of the lack 
of relativistic covariance of formulae (4.10). However, this lack is only apparent 
since we may replace the particular time instant & by a particular time-like hyper- 
surface o, which is a relativistic concept. 
— — 

з As j, depends on 2) , this > properly speaking, not a solution but a transition to an integral 
e quation (which may be solved by iteration). 
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Let us consider the operator of the energy density (4.1) at the time instant tg 
and introduce the solution A, in the explicit form (4.10). Taking account of (1.5'^), 


we have 
pas р IT : 
Hit) = H,, + HO + un A р + grad AN ° grad Ap + gan HG 
1 EE (4.12) 
LE grad 469. grad A® — Am], — 460 j, + т (A + A) k, 


where all the field quantities are taken at the time instant tọ. The operator НО), 
is the same as (1.5’’) with A, replaced by сел» апа represents the energy density 
of the photon field. 

Going over from the energy density to the total energy, we may perform some 
integrations by parts, e.g., 


| x grad 4% grad A = —f dx AD A 4? = f dx ue (4.13) 


whereby some of the terms cancel and we are left with 


A р ИТУ : 1 f 
H(t) = Hy, + H9, + [es Е OAP EAO ADEA), 


1 
)4.14( | 1 )40 + سی 2+ 


Although the total energy of an isolated system of fields is constant in time, we 
have written Н (tj) in order to remember that all the quantities in (4.14) are taken 
at the time fp. | 

Let us denote 


ie = 4 Р 
IE № 2 (—y, grad y, + grad р," Pa) | (4.15) 


From (4.3) it is seen that 
j = 9 + в (4.16) 


and (4.14) may be written also in the form 


на) = Hein, + НӨ), + [| аз (ADAP + AOA 


1 1 
e (st) ;(0) — 400) 
о 4 ii Ir 2I (4.17) 
As a matter of fact, this expression is extremely complicated since, in our case, 
the current (and also the potentials Аб) depend on A „Which may be iterated again 
and again so that the energy represents a power series in the coupling constant. 


\ 
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By limiting ourselves to second order terms, and taking account of (4.11), we are 
left with? 


1 д (0) у. (0) уу 
H (to) = ШИ un H, m. E I 94, (ғ) 9h 6) PTO eru 
Arc 9% 91 | PPL үл | 


д ðt | x gt! | 


t 
t3 Е DES Ex. eT | РЕ mu 


1 с 1 1 
, rs (0) „ус 0 , , ГА 
rE [es ТА АҒУЫ (x^) j (x E me = = fe AM(26')kO(x') (4.18) 


The first two terms are of the order zero in the coupling constant. The third term 
is of the first order and should be treated by means of the (time independent) 
perturbation method. The remaining terms are all of the order e?. We distinguish 
among them two types of terms: (2) the cross terms (in which the current of a field 
is multiplied by the current of the other field); they describe the direct interaction 
between the two fields, and (fi) the remaining terms. If only one particle of each 
type is present, then the remaining terms are of a self-energy type (an interaction 
of the particle with itself). We shall not deal in this paper with the self-energy 
terms but we shall confine our attention to the cross terms. These terms may be 
picked up from (4.18) and are (4.19") 


Ат d 1 [м — x’ Kies 


T‏ ا 
ME (4n)? |x’ — ж” | | x | ,‏ 
pa AR M («”)‏ 2 


as er Fer ЕДІ 377 р 
T (дт)? | x’ | ті о,“ 


where the indices 1,2 apply to the two charged fields with masses m,, m, and 
charges e}, еҙ. Up to terms quadratic in the coupling constant the Hamiltonian is 


H =H; +A+B+C+D+H+ THOU (4.20) 


(4.19”) 


where H®, may be simply omitted if we consider a state without photons, H 


4 We denote O — (110), 1, 40), А00) = 0 where & is obtained from & by replacing A, 
by AP in (4.3) x’ stands for the four variables x’, t. 
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МЛ. ^“^ __ тз с GG nennen 
is the term of the order е together with all the self-energy terms (of the order е?) 
and may be considered as a perturbation. The remaining terms are the kinetic 
energy and the cross terms (4.19). These terms may be treated in the configuration 
space. Consider the representative < | хі; xı) of the state vector in the subspace 
of one particle of each type (x, is the coordinate of the electron, and x4 of the 
proton) and write the equation 


(|Ң, + 4 +B + C+ Dl х;у = W |х; ху» (4.21) 


The rules (2.22) —(2.25) applied to the two types of variables x, and x, indepen- 
dently enable us to transcribe (4.21) into the configuration space. With the aid 
of (2.22), we get immediately 


д 
a (misa (4.22) 


ЫН 


2 
1 
LH. ; = Яз: 
<| kin. |х х1) = E. 
Rule (2.24) applied to the term A yields the Coulomb interaction 
e 
(| 4 | хи; x = enin (4.23) 


where г denotes | x, — х, |. For the evaluation of B, C and D the application of 
rules (2.25) is required. We get 


< В yas ете» 1 ОО а: di 9 l д 
хх.» = — 
т is түт, [г 9% Quia ont т] д& 


1/19 16% 1-459 (4.24) 
+5 (= uta er E aoe 


This term describes a stationary interaction of the two currents and vanishes in 
the limit of an infinitely heavy proton (mg œ), which is obvious since such a pro- 
ton would rest and produce no current at all. Taking account of the identity 


1 
A =—476° 
à To (к) (4.25) 
the last term may be written also in the form of a ”contact interaction” 
1 NN 
6) 
4; EE (r), (4.26) 


which vanishes unless the coordinates of the two particles x, and x, coincide. 
In order to compute the term C we invoke the continuity equation (1. 4) 
whereby (4.21) 


1 
С==— ке ae d v | n 
a olf er 
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The application of rules (2.25) yields 
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We notice that the last three terms are identical with the last three terms of (4.24). 
The integral over d?x in the first term may be evaluated in an elementary way 


Ie 1 ) 5 1 
ae) a) 62» 


where 
1 (x —^5) (хі — x) 
T; = e б,, ТЕ қта қ лын (4.30) 
апа 


r= | xı — x; | = (х —х) (x — x). (4.31) 


Both terms B + C yield together 
ете, b 9 9 12109 9 


1 
(|B C|x >= — 


4т mm; 


h ك‎ д + ОЕ д + 2x8 (r) 
9x т) әж эт осы Е CRATE 
The evaluation of the term D given by (4.19) is rather lengthy. To get rid of the 
time derivatives of the field quantities in D we take account of the field equattions 


(1.7) in a simplified form (without interaction terms) 


о O т 0 dak 


(4.32) 


Ay, = Wa | (4.33) 


2m 


as we are interested only in term not higher than quadratic in the coupling con- 
stant. A lengthy but elementary computation yields 
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5. The transition to relative coordinates 


We introduce the relative coordinates, and the coordinates of the mass centre 


тух + Myx 
x= ухх, X ET MOL Ui (5.1) 


m, + тэ 


and look for a so ution corresponding to the mass centre at rest. Thus, the repre- 
sentative < | xı; х1> should be independent of X. In this case equation (4.21) 
simplifies and may be written as a Schródinger equation in the relative coordinates 


1/1 1 ее; 1 1 ее. luo 
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= Wy (x) (5.2) 


where T; is given by (4.30). The first term describes the kinetic energy with the 
reduced mass constant as usually, the second term is the usual Coulomb inter- 
action (unipole), while the remaining terms describe corrections due to the rela- 
tive motion of the two particles. Some of the corrections are of a dipole, some of 
a quadrupole character. We notice also the appearance of contact terms. These 
corrections are all extremly small in the case of a hydrogen atom since they are 
roughly of the order A,A,/a? or 724 а respectively, where A = 1/m is the Com- 


pton wave length, while a is the atomic diameter. However, similar terms іп а corres | 
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ponding theory of the deuteron (coupling, e.g., by a vector meson field) may 
be of great importance owing to the fact that the radius of the deuteron is roughly 
of the same order of magnitude as the Compton wave length of the nucleon. 


КРАТКОЕ СОДЕРЖАНИЕ 


И. Райски, Енергия связанных состояний в квантовой теории полей. 


Из основных правил теории полей выведено уравнение собственных 
значений энергии системы двух тел. Для примера рассмотрено простой модель 
не релятивистских, не имеющих спина, заряженных полей. Так как про- 
блема соответственного значения энергии соотвествует установленному мо- 
менту времени, можно избежать рассматривания опоздания, а временное 
электромагнитное поле можно изобразить в виде суммы двух выражений: 
первое — в виде свободных волн, другое — в виде стационарного решения. 
Оператор энергии поля расделяется на ряд частей: кинетическую энергию 
заряженных полей, энергию фотонов, энергию взаимодействия между Фото- 
нами и заряженными полями, разные выражения собственной энегрии (взаи- 
модействия заряженного поля с электромагнитным полем, вызванным тем же 
заряженным полем), выражения смешанные (из взаимодействия заряжен- 
ного поля с стационарным, электромагнитным полем, вызванным иными за- 
ряженными полями). Выражения смешанные могут быть записаны в конфи- 
турационном пространстве и составляют источник (кроме куломбовского вза- 
имодействия и выражений Брейта) некоторых добавочных поправок. Урав- 
нения собственных значений энергии, т. e. обобіценное, независимо от времени, 
уравнение Шредингера, написано в относительных координатах двух частиц. 
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THE INFLUENCE OF HINDERED ROTATION ON THE 
SCATTERING OF SLOW NEUTRONS BY BOUND PROTONS* 


By WEODZIMIERZ Koros 


Institute of Theoretical Physics University of Warsaw, Warsaw 


(received September 18, 1953 ) 


The cross section of the proton bound in a molecule which has an additional degree 
of torsional freedom is calculated. It is shown that the cross section of the proton for 
the collisions which are elastic with respect to the hindered internal rotation but exciting 
the rotation of the whole molecule about the axis of the hindered rotation, is smaller 
than that of a proton bound in a rigid molecule. It follows from the theory that the low- 
ering of the cross section depends on the energy of the incident neutrons and on the 
thermal energy of the scattering molecules. 


l. Introduction 


So long as the neutron energy is large compared to the oscillation energy of 
the proton in the molecule we may consider the proton as free, and its cross section 
will be nearly independent of the energy of the incident neutrons. This condition 
is not fulfiled by thermal neutrons and it was Fermi (1936) who first calculated 
the cross section of the proton for the scattering of slow neutrons taking into 
account the influence of the oscillations of the proton in the molecule. Fermi 
showed that for sufficiently slow neutrons the scattering cross section may be 
calculated assuming the interaction potential between neutron and proton to be 
of the form of a delta-function and using the first Born approximation. More 
detailed calculations in this approximation allowing for lack of symmetry in the 
binding, thermal motion of the scattering molecules, etc., һауе been carried out 
by Bethe (1937) and Arley (1937). 

More exact treatments of the same problem by an integral equation method 
given by Breit (1947), and by the variational method of Lippmann and Schwinger 
(1950) given by Breit and Zilsel (1947), Zilsel, Darling, and Breit (1947), and 
Lippman (1950) yielded only very small corrections of order 0,3 per cent, or 
less, to Fermi's result, thus showing; at least for some special cases and very slow 
neutrons, that Fermi's approximation is a very accurate one. 

Sachs and Teller (1941) considered the case when, the neutron energy being 
small compared to the quanta of molecular vibrations is nevertheless great com- 
pared to the energy differences between the rotational levels of the scattering 
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molecule. Replacing the proton by a freely moving mass-point whose mass is a ten- 
sor depending on the mass and structure of the molecule, they have been able 
to calculate the cross section of the proton bound in a molecule. The theory for 
the case when the thermal energy of the scatter- 
ing gaseous molecules approaches zero, leads 
to a very simple result: 


1 
с = б (Mla): (1) 
where o, is the cross section of a proton rigid- 
ly bound to an infinite mass, and the џ;'ѕ are 
connected with the mass tensor of the proton. 
For some scatterers used in the liquid state 
the proton cross section calculated from for- 
mula (1) is in very good agreement with experi- 
mental data, as may be seen from Table 1. 
We see, however, from Table 1 that this 
agreement is only valid for the simplest molecules. 
In the case of methyl alcohol (CH,OH) the 
Fig. 1. Model of the methyl alcohol experimental value of the cross section found by 
molecule Janik (1953) differs considerably from that given 
by formula (1). 
It is natural to expect that the internal hindered rotation which can 


take place in the molecule of methyl alcohol about the axis C—O (Fig. I) is the 
cause of this difference. 

In Section 3 it is shown that the appearance of internal rotation lowers the 
proton cross section for the collisions in which the rotation of the whole molecule 


about the axis of the hindered rotation is excited. In the last section we apply our 
result to the molecule of methyl alcohol. 


Table 1 
Scatterer б per molecule с per molecule 
calcul. from (1) Са References 
48,6 . 10-24 cm? 48 . 10-24 cm? Carrol (1941) 


H,O 78,5 . 10-24 cm? 81,9. 10-24 cm2 Janik (1952) 
CH,J 195 . 10-24 cm? 195 . 10-24 cm? Janik (1953) 
CH,0H 231 .10-24 cm? 180 .10-24 cm? Janik (1953) 


2. Fermi's approximation 


~In Fermi’s approximation the differential cross section of a proton for a col- 
lision in which the proton is excited from its m; -th state to the m,-th state is given by 


ке 
т К.т 


l 
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where U, ‚(r,) and О, 749 г) are the eigenfunctions of the proton in the initial and 
final state respectively, 1 m, and m, denote each a complete set of quantum numbers 


for this states, r, and r, are the coordinates of the proton and neutron respectively, 
oo is the cross section for a free proton, u’ = М ot ММ, һы + М, is the reduced 


m 
mass for the relative motion of neutron and molecule, 


27 \? 20.5. 2u 2u'E 
k2 — E i > рТ; деті. em e Pr и Pg 
i ( 4; ) h? vi h2 E Ey — Em) = h^ 


EM " (3 
kif = a kr, и = М 9 


E; and Е; denote the energy of the neutron before and after the collision respec- 
tively, E, and E, аге the energy of the proton in the initial and final states res- 


pectively. The difference between the neutron and proton masses is neglected. 
Having integrated over the coordinates of the neutron and chosen a coordinate 


system with the z-axis in the direction of kj, we get 


Клио NA ج‎ 2 
трт T EH | exp{tkz,} Un (ть) U mr) Wp dQ (4) 
where for abbreviation we have written k = | kgl. 


For very small neutron energies, when all the scattering molecules are in 
_ their ground states, only elastic collisions are possible. In this case, replacing the 
exponential function in (4) by unity (k — 0) and integrating over the full solid 
angle, we get for the total cross section 

lim 060 = Au?o, (5) 
Е; >» 0 
as expected. | 

Using formula (4), we may calculate the cross section of a bound proton. 

But in the case when we want to consider the influence of the rotation of the scatter- 
ing molecules the method becomes impracticable. The differences between the 
rotational levels are small compared to the thermal energy of the molecules for 
the temperatures commonly used in practice. We must, therefore, take into 
consideration the statistical distribution of the scattering molecules over the 
rotational states, and then calculate the cross sections for all possible transitions. 
The cross section which can be compared with experiment is the sum of all to- 
gether. This requires a calculation of a very great number of contributions and 
so the method becomes in general impracticable. But nevertheless we shall show 
that this method enables us to get in a simple way some conclusions concerning 
the influence of the hindered rotation on the proton cross section. 


3. The influence of hindered rotation 


The molecule will be taken to be a symmetric rotator which in addition to 
the usual three degrees of rotational freedom has about its axis of symmetry also 
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E 


a degree of torsional freedom between two of its principal parts. The motion of such 
a molecule is best described by means of the Eulerean angles 9, Y, фі, ys, where à 
denotes the angle between the space fixed axis z and the axis of symmetry z' of 
the molecule, y is the angle between the line c* nodes and the x axis, and 9, 


and фо are the angles between the line of nodes and the х — axis fixed in the 
lower part of the molecule and the x’ axis fixed in the upper part of the molecule 
respectively. 
We assume the hindering potential to be of the form 
H 1 
Va m (1 — cos n 9^), (6) 


a 


where His the height of the barrier, y’ is the angle of rotation of one part of the 
molecule with respect to the other, and 7 is the number of minima of the potential 
energy curve as g’ increases from zero to 2m. The potential function given above 
is one of the simplest having the desired flexibility and necessary symmetry prop- 
erties. 

By introducing new variables defined by 


ф' = 91 — Ф», 9 = В.Ф. + Boa: (7) 
where 


ß ái (125172) 8 
= — == " 

1 A, 4, | ( ) 
and A, is the moment of inertia of the i-th part of the molecule about the axis 
of the internal rotation, the corresponding Schrédinger equation has been found 
(Nielsen 1932) to be separable into two differential equations, one being the quan- 
tum mechanical equation for an ordinary symmetric rotator, the other having 
the form of a Mathieu equation. 

This leads to eigenfunctions of the form 


1 
Eg; №, Ф, 9^) iT) Von bs (9, №, Ф) Mos (o^). (9) 


where о is the quantum number of the hindered rotation, m denotes a set of quan- 
tum numbers (r, z' and j) describing the rotation of a rigid molecule, U,,(8,y,9) 
is the normalized wave function of a symmetric rotator given by Reiche and Rade- 
macher (1926), and Kronig and Rabi (1927). 

U, (9, y, v) = C', Oder tr (10) 


where 
stb gs 1 
6 „(9) = 22 (L— t)? F(m;t) t =5 1 — c4), лз ҮЙ) 


Е(т; 1) =F (Ê. )ئ‎ is а hypergeometrical function with arguments depend- 
ing on the three rotational quantum numbers 7,7’, j denoted by m, С’, takes 
care of the normalization, d and s are defined by 


4 = | +1'|,s=|7—1'| (12) 


— = 
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M ,m(p') in (9) is the Mathieu function which according to Floquet’s theorem 
(Whittaker and Watson 1950) may be written in the form 


+ co 
Mes OD se Er y a eer аз) 
n--oo 
Letting фу and g, increase independently by integral multiples of 27 and 
using the condition of single-valuedness of the wave function У 00,9, v. $’), 
we can calculate the characteristic exponential of the Mathieu function: 


с = —t f. (14) 
The coefficients a, of the expansion of the periodic part of the Mathieu func- 
tion in a Fourier series depend on о and on the characteristic exponential с, i.e., 
by relation (14), опт, the quantum number of the rotation of the whole molecule 
about the axis of the hindered rotation. 
The z-component of the radius vector of the proton in terms of the introduced 
angles is given by 
z,, = €, cos 0 + a4 sin 9 sin (p + Bap’) (15) 
for the proton bound with that part of the molecule which has the moment of 
inertia A, and whose rotation is described by the angle ф, and 
2,, = а соз Ü + аҙ sin 9 sin (y — fg’) (16) 
for the proton bound with the other part of the molecule. The a; and c; in (15) 
and (16) are connected with the geometrical dimensions of the molecule. For 


8 =y =p —9' =0 we have x, —aj; у, =0; 2, = с. (For the methyl 


Pi 
alcohol molecule see Fig. 1). 
For the calculation of the differential cross section we substitute (9) and (15) 


or (16) in (4) and get 


da Ж) 


малта” m)? Готама cos 9 + agsin 9 sin (y — Р19')]} Um, (v.g) 
a :7(277 


Г d2. oy 


x U,, (9. ¥, Ф) Mom; (9) Мет; (P) dv 


The product of the Mathieu functions may be written in the form 


+оо 2590 сад 
M (Ф)М, | (ф') = eei — of) g' У Жаса - е— т: — yh pog : (18) 
ا‎ eimi Je 


п= —oo 


where 


oo 

ь = V aen? ger? (19) 

n m m-n 
m= —oo 


Substituting (18) in (17), using (10), and changing the integration variables: 


А р — Pp’ = 9 
s (20) 
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we get 
л 2л 2n 2n 


ANT ES Ce b, f expfik(c,cos® +a, sin sind! 
my ofymioi rer)? Со, т; » fff pt (со 2 ) 
n=-© 0000 
d; жар Si tse 
+ i(t,—t,)d + Кт; --т/)у% ар (Тар ы 
2 


x F(mj; t)F(m;; де?” sin 0 düdyd9d9' | dQ. (21) 


Ж 
do en 


Integration over 9” gives at the right hand side as the only non-vanishing 
term that for n = 0, and we see at once that the resulting expression is 52 multi- 
plied by the differential cross section of a rigid molecule, i.e., the cross section 
which may be obtained from (17) by putting ф’ = 0. 

Thus, we may write our result in the form 


do,, 2. ас т. (22) 
where do,, mi is the differential cross section of a proton bound to a rigid rotating 
molecule. 

For the total cross section integration over the full solid angle gives an ana- 


logous relation: 
иней Е О (23) 
The formula (23) gives the cross section of a proton bound to that part of the 
molecule whose rotation is described by the angle gg. For the proton bound to 
the other part of the molecule we must substitute in (4) (15) instead of (16). But 
the Mathieu function may be written in the form 


+ со + со 
М mlp’) МЕ ea thie" 3 а ЛУ = eit hae’ 3 a, e" ae" (24) 
n- —oco п-ш--со 
and the product 
Too 
* Қт;-ту/у)р p ino? 
Сы Toe (25) 
n= —оо 
where 
+00 
= (0;т;) gg (erty) 
2 deitate etr (26) 


п=— со 


Now, using (15) and (25), we get. from (4) also the result (23) for the total 
cross section. In this case b, is given in terms of the same a,, which are only TS 
with other indices. 


From the normalization of the wave function Pino (9), we have 


Too 
dary =1 (27) 


n=—00 
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Too 
У @ep)t=1 (27) 
п=— со 
which, according to Cauchy’s inequality, gives 
+оо 
b, = У асас 
п---со 


From (27) and (28), we see at once that the cross section for collisions in 
which neither the hindered rotation nor the rotation of the whole molecule about 
the axis of the hindered rotation is excited will be the same as for a rigid molecule 
(by = 1), and for collisions in which the latter is excited, we shall have a lowering 
of the cross section compared to that of rigid molecule. 


4. Applications 


In the preceding section we have been led to the conclusion that the proton 
cross section for inelastic collisions with respect to the rotation of the whole mole- 
cule about the axis of the hindered rotation is smaller than that of a rigid mole- 
cule. Now we want to estimate the magnitude of this lowering for the molecule 
of methyl alcohol. The molecule will be taken to be a symmetric rotator. It will 
be assumed that the axis C—O is the axis of symmetry. The slight asymmetry of 
the actual methyl alcohol molecule is due to the relatively light hydrogen atom 
in the OH bar, and Borden and Barker (1939) have shown that the two moments 
of inertia of the methyl alcohol molecule about the x and y axis perpendicular 
to the axis C—O differ only by 4 per cent from their average value. 

Koehler and Dennison (1940) have calculated for methyl alcohol the coeffi- 
cients а, of the Fourier expansion of the periodic part of the Mathieu function 
for a definite value of the barrier height in the hindering potential (6). Using their 
results, we may easy obtain the values of b? for various transitions. Some of them 


For the proton 
bound in OH 


are given in Table 2. 


Table 2 


See 
Fox [оз [ов [ur [oro 


For the proton bound in the CH, group. 


Transition 
<-> 
т; т, 


The arrows indicating the transitions have been written in both directions, 
1.е., for ordinary inelastic and for so called ”superelastic” collisions (Е, < Ey). 
This is allowed because bọ depends only on the initial and final states of our system, 
and it may be shown that there is a kind of symmetry between these states, and 
that their interchange does not influence the values of bọ. | 

From Table 2 we see that for the collisions (elastic with respect to the hin- 
dered rotation) associated with transitions consistent with the optical selection 
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rules, we have only a very small influence of the hindered rotation on the proton 
cross section. The experimental data, however, show this influence to be of the 
order of 22 per cent. On the other hand b? is smaller, i.e., the influence of the 
hindered rotation greater, for higher "multipole" transitions, and, in analogy 
with optics, formula (17) may be regarded as collecting all "multipole" transi- 
tions, which may be obtained explicitly by expanding the exponential function 
under the integral sign in a power series. A comparison of the neutron energy 
with the energy differences between the rotational levels of the methyl alcohol 
molecule shows that, for thermal neutrons, the former, is sufficient to excite the 
high "multipole" transitions associated with a large change of the quantum 
number т, e.g. 0-8, 1—9, etc. It is also justified to expect a relatively great prob- 
ability for such transitions because the neutron wave length is of the same order 
of magnitude as the molecular dimensions. 

The theory of Sachs and Teller (1941) gives a relationship between the proton 
cross section and the neutron energy, and formula (1), for gases, holds only in 
the limiting case when the energy of the incident neutrons is much greater than 
the thermal energy of the scattering molecules. We may expect, however, the 
effect to be small for liquids, which follows from the good agreement of formula (1) 
with experiment. 

However, from Tab. 2, we see that 52 approaches unity with decreasing | 


so we may expect that the cross section of a proton bound in a molecule such as 
CH,OH will depend considerably on the neutron energy and on the thermal 
energy of the scattering molecules (the latter owing to "superelastic" collisions) 
The lowering of the cross section caused by the hindered rotation will be smaller 
when both the neutron energy and the thermal energy of the scattering molecules 
decrease, and when they approach zero we obtain the same cross section as for 
a rigid rotating molecule. 

From the experiments published up to now, it is not evident if the expected 
effect really exists. I am, therefore, very indebted to Dr J. A. Janik for his promise 
to carry out measurements which will allow to settle the question. 

I should like to express my sincere thanks to Professor L. Infeld for his kin 
interest in this work, and to Dr J. Werle for many valuable and helpful discussions 
in the course of the calculations. 


Note added in proof. 


The coefficients bj given in Table 2, are calculated assuming the height 
of she potential barrier to be H = 760 cm-!. However, E. V. Ivash and 
D. M. Dennison (J. Chem. Phys. 21, 1804 (1953)) using recently published 
spectroscopical data on the microwave spectrum of methyl alcohol have found 
the height of the potential barrier to be only 374,8cm-!. The coefficients bà 
caleulated with this new value of the potential barrier are smaller than those 
given in Table 2. E. g., for the proton bound in the CH, group we shau hawe now 


Scattering of Slow Neutrons by Bound Protons 15 


HT 0=->1 0«——4 ОТ | 010 | 


It 


bs 0,98 | 0,88 0,75 0,58 | 


КРАТКОЕ СОДЕРЖАНИЕ 


В. Колос, Влияние задержанноо вращения на рассеяние медленных нейтро- 
нов связанными протонами. 


Вычислено эффективное сечение для рассеяния медленных нейтронов, 
протонами, связанными в молекуле, обладающей прибавочной степенью тор- 
сионной свободы. 

Доказано, что эффективное сечение протона на столкновения, которые 
являются упругими относительно заторможенного внутреннего вращения, но 
неупругими относительно вращения целой молекулы вокруг оси внутреннего 
вращения, — меньше чем эффективное сечение протона связанного в жесткой 
молекиле. Из теории следует, что эффективное сечение зависит от энергии 
падающих нейтронов и термической энергии рассеиваюших молекул. 
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ON THE MASS SPECTRUM OF ELEMENTARY PARTICLES 


By Jerzy Ravski 


Nicholas Copernicus University Toruń 
(received November 18, 1953) 


Generalized field equations (Schródinger-Gordon and Dirac equations) together 
with reciprocal supplementary conditions postulated within the framework of the non- 
local (bilocal) field theory are shown to lead to mass quantization. Families of particles 
may be classified with respect to (integer or half integer) values of the spin and of the 
isotopic spin. The mass spectra of the nucleon and pion families are given. The masses 
consist of two parts: a proper mass and a field mass.The existence of a field mass is re- 
sponsible for a fine structure of the mass levels. 


l. Introduction 


In view of the experimental evidence of existence of families of elementary 
particles the problem of a theoretical explanation.of their mass spectra becomes 
a matter of urgency. We present here an attempt to forinulate a field theory in- 
cluding mass quantization based on the idea of bilocalizability. 

Let us begin with a general remark. Two classes of quite different characters 
may be distinguished among physical quantities appearing in the field theory. 
The position and displacement operators x,, d,, the spin variables (or indices) a, 
and the isotopic spin variables (indices) f belong to one of the two classes. On the 
other hand, the field quantities y and all the operators defined as functionals of y 
belong to the other class. The quantities belonging to the first class are universal 
while those belonging to the second class are specific. Indeed, the operators x ,, 
denoting position (and time), and the operators id, denoting the momentum (and 
energy) are universal since they do not specify what kind of particles assumes 
the positions and momenta given by the eigenvalues of these operators. On the 
other hand, the momentum operator constructed with the aid of a field quantity y 
refers to that kind of particles which are described just by ф and se it is a specific 
quantity. x,, d,, a, В together form a field of a geometry in a broad sense of this 
word. Owing to the inclusion of d, it is a dynamical geometry. 

From the point of view of the above classification a defective character of the 
traditional local field theory is evident: x,, a, and f are arguments of the specific 
field quantities y, while in a complete field theory the field quantities are expected 
to depend upon all the universal quantities HI 

ү = р(х, 4„ а, В) = Фк,» d.) = р(х, d) (1.1) 


(77) 
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where 
[x,, 4] =д,,. (1.2) 
In this way the "non-local" fields are viewed as an almost unavoidable generaliza- 
tion of the traditional "local" field theory. 

In consequence of (1.1) the field quantities become operators in a twofold 
sense: as quantities subjected to field quantization and as functions of the opera- 
tors x, and d,. In the x-representation the operators x, and d, as well as the field 
quantities y become matrices: 

2 x , / 9 2 
«x xx Se (а a ed а = ed о (1.3) 
и Hu ۳ ay 
и 
y = «x | p(x, d) |x" > (1.4) 
satisfying the matrix algebra. Sometimes it is convenient to regard y as function 


of Yukawa (1950) variables 
y = y (х, т) with х= 5 (x, + x,) pig es (x, — x,)- (1.5) 


From (1.4) and (1.5) it is seen that the name "bilocal" fits better than "non-local" 
field theory. (The name ”bilocal” field theory was introduced by Bagge 1953). 

There exist two independent linear forms in y and d,, namely their commu- 
tator and anticommutator. In the x-representation we have 


Iyla;r) 
erts SEEN ew e E 
x А = | — ed ў == 5 
ENT 9х 9х Е Е ylar) 450, 

dr, 


There is a close analogy between the two correspondences: classical <-> 


quantum theory, and local <-> bilocal theory. In the local "limit" w(x, г) becomes 
independent of r, so that 
9v() 
<<“ [rar 1% OC), (1.7) 
0. 


Thus, the anticommutator [y, d,], has no analogue in the local field theory and 
may be used for introducing new features into the field formalism. 


2. A generalization of the wave equations 


In the interaction-free case every local field quantity satisfies the Schrédinger- 
Gordon-Klein equation so that the d’Alembert operator 


ә 9 9? 
re ee С хр 2.1 
ie: Ix, ex. че ді (2.1) 


appearing in it acquires a fundamental significance. A question of great impor- 
tance is what is the proper generalization of the d’Alembert operator for the bilocal 


case. There are two such independent operators symmetrical with respect to. x 
nad i 
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9 9 22 22 
MATT 2 
des CE 2 9472 em 
and any linear combination of them 
4С ARG 2C uM да 
Diese trie ь iid 
u и и 
with 
By putting 
1—C 26 
С = —— ber к= 
= ‚= 5 (2.5) 


the operator (2.3) may be written as 


ка tar) (2 ==) 


with an arbitrary dimensionless constant C. Going over to Yukawa variables (1.5) 
the operators (2.2) are respectively 


9? 9? 9? 9? 
xi = oF and Oa v: 22 , (2.2) 
while in the general case 
9? 9? 
00 = E +4 ЕТ 2 (2.6 ) 


With the aid of the generalized d’Alembert operator we may postulate a generalized 
wave equation 
ПО y = 0. (2.7) 
This wave equation for an interaction free field has been written down without 
92 
‘the usual mass term ..? since the operator 4C gx may be just considered as a sub- 
A 


stitute for the mass parameter. Equation (2.7) may be written also in the operator 


form 


Пу, dl- , а,]- TE C[[v. а,] + ? d l+ =0 (2.8) 


which follows easily from (2.6) and (1.6). Three cases have to be distinguished 
C>0,C<0 and С = 0. In the case С = 0 we fall back into the local field theory 
of particles with vanishing rest mass while in the case C + 0 the eigenvalues of 

22 


the operator 4C Ore 


are expected to determine the mass spectrum of families of 


articles. [Jw 
р Similarly as in the local theory the wave equation (2.7) may be linearized. 


i i i ation 
Let us introduce two sets of matrices a,, f, and write the first order equ 


aly, d,]- = ia Bly 4,] + (2.9) 
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or 
9 9 хк. NM С А "hay | x" (297 
„(э + ge (< =, (g7 ¬7 (>< @ 
p и 5 1 
Or 
а, — P(X г) =21ар, 7 фиг : 
А x, lad Or, 


where a is a dimensionless constant. Let us assume e.g. that a, and В, are two 
sets of matrices setisfying the anticommutation relations 


la,» а] + = Ж Bila = 26, , [а„, ЖЕ == 0; (2.10) 
Multiplying (2.9’’) by a, ox We have 
9? ә 29 д д 


ont vy (x, г) = 2aia,B, 2 or, p(x, т) = -2aip, EA Е y (xr) 


ә 9 2 
= — 4a? , E» dr y (x, г) = —4a? 22 (zT). (2.11) 


From this it is seen that the generalized Dirac equation (2.9) may be carried over 
into the generalized Schrédinger-Gordon-Klein equation (2.8) with a positive 
C = —а?. The case of a negative C cannot be excluded either. It is immediately 
seen that a reversal of sign of C is obtained by assuming that the two sets of matrices 
anticommute 


[2,, Ble = 0, 2.10) 


Now, there arises the question of the physical meaning of the degrees of 
freedom connected with the matrices a, and ß,. The a,’s are obviously connected 
with the spin of the particles in exactly the same way as in the local theory. Assum- 
ing the standpoint represented by Pais (1953) who pointed out the necessity of 
an intrinsic incorporation of the isotopic spin into the framework of the field 
formalism, we believe that the Bis are connected. with the isotopic spin. 
There exist four possibilities: (1) half integer spin and isotopic spin, (2) 
integer spin and isotopic spin, (3) half integer spin and integer isotopic spin, (4) 
integer spin and half integer isotopic spin. The first possibility is realized by 
nucleons, the second by pions, the third probably by leptons, the fourth is yet 
unknown. 

(1) We assume that the nucleonic field is described by the generalized Dirac 
equation. (2.9), where a, and В, denote Dirac matrices satisfying the relations 
(2.10). In this case the Dirac equation may be carried over into the wave equation 


(2.6) with a postive C — аў, ‚ where the suffix + has been attached in order 


А 5 . sgh 
to remember that the constant а, applies to the case of the isotopic spin = 
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(2)We assume that pions are described by the second order wave equation 
(2.7) or (2.8) with a (positive) constant C == +a, where the suffix 1 has been 
attached to remember that the constant a, appliesto the case of the isotopic spin 1 
Now, the question arises as to the relation between the constants ay, and a, which 
cannot be accidental. The simplest possibility would be dy = а but the compar- 
ison of the mass spectra of the two families would lead in this case to a complete 
disagreement with experiment. The next plausible assumption is that their ratio 
is connected with the dimensionality 4 or 10 of the irreducible representations 
describing spin 1 and 1, or more precisely with the number of independent һурег- 
complex numbers 4? and 10? respectively. We assume therefore 


a 
Аға, = 10а: 50 that а зе а = س‎ · (2.12) 


2 


(3), (4) The cases of disconnected spins, i.e. of a half integer spin and ап 
integer isotopic spin (or vice versa) make plausible the assumption that both 
sets of matrices а, and В, commute with each other and that the а,5 аге Пігас 
matrices while the ß’s аге Duffin-Kemmer matrices (or vice versa). The correspond- 
ing second order wave equation (2.7) or (2.8) assumes a positive constant C > 0. 

Besides these cases there is a local case C -- 0 having no connection with 
the isotopic spin whatever, and realized by photons and probably by gravitons. 
Hovewer, it is not excluded that these particles are contained in the case of integer 
spin and half integer isotopic spin. 


3. The supplementary conditions 


Examples of mass operators of a similar type as that given by (2.9) 
2 (3.1) 
М = +2aß, ar, .1) 


have been investigated recently by Yukawa (1953). However, the eigenvalue 
problem for such (fourdimensional) operators meets with serious difficulties 
connected with the non-definite character of the metric form. Indeed, dealing 
with the eigenvalue problem in exactly the same way as one used to do in the older 
quantum mechanics with three-dimensional operators (i.e. operators which are 
scalar products of threevectors) Yukawa found expressions contradicting the 
Lorentz invariance. The reason for this difficulty lies in the fact that we cannot 
treat the fourth component of a four-vector on the same footing as the three space 
components. Above all we cannot assume boundary conditions (e.g. vanishing 
of the wave function for ry — = oo). Time-like components cannot be used for 
defining discrete eigenvalues. Thus, any formalism involving mass operators in 
the form of scalar products of four-vectors is too wide and has to be limited so as 
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to restrict them effectively to space-like four-vectors. To this end we need some 
supplementary conditions}. 
Most probably the character of such supplementary conditions is universal. 
In other words, the supplementary conditions should be the same regardless of 
the character of the field quantities. Thus, the matrices a,, f, are not expected 
to be involved into the supplementary conditions. Moreover, universal conditions 
(restricting the effective range of variability of r,) being a matter of geometry are 
supposed to possess covariance properties characteristic (intrinsic) for the geometry 
itself. Now, our generalized geometry involving x, and d, is characterized not 
only by the Lorentz invariance but also by an invariance against the transformation 
of reciprocity (Born 1949) 
Atx, => iAd,, iAd, > Mix,» (3.2) 
where the constant A has been introduced for dimensional reasons. (With c = h = 1 
the constant A is of dimension of length). It seems plausible to require that the 
supplementary conditions be invariant against (3.2). Such conditions introduced 
previously by the author (1951,a, b) are? 


[[v. а. ЖЕ = |у, Xl- d,]- = 0, (3.3) 
Ally, d,]-, d,]- — 27Ір, х„]-, x,] = 0, (3.4) 
With the notation 
sis 1 1 NR 
UE Ys (ЛЕ х, | 1 Àd) (3.5) 
(3.3) and (3.4) may be written in an equivalent but more elegant form? 
(ly. at]-, af]. = lly, az]. az]. = 0. (3.6) 


It should be remarked that the supplementary conditions (3.4) or (3.6) are recip- 
rocal only formally since the field quantities y are not reciprocal themselves 
(as the wave equations (2.8) or (2.9) are not reciprocal). This seems all right since y, 
being a specific quantity, must possess properties reaching beyond pure geometry. 
Only the constraints determined by the ln conditions but not the field 
quantities need to be reciprocal. 

With (3.2) a constant of dimension of a length A has been introduced into 
the formalism. It is not obvious whether A is a universal lenght or whether it is 
a "rate of reciprocity” with a different numerical value for each type of field. 
At the present stage of the investigation we have not been faced with the necessity 
of assuming different 4's, so that we may regard it as a universal length. 


1 Supplementary conditions exist also for the local field theory, namely 
[v, x)= 0, WHS 1,...,4 
satisfied by any field quantify independent of d. — 
2 Yukawa (1950) assumed [[y, x,] , x,]_ — 4% instead of our (3,4). This condition is simpler 
but not self-reciprocal. 
3 A more general form with reciprocal invariance may be obtained by replacing zero to the 


right hand side of (3.4) or (3.6) by Cy. We have put the dimensionless constant C equal to zero to 
avoid new, not absolutely necessary constants. 
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In terms of Yukawa variables the supplementary conditions are 


д 595 
0 a ci uen co (3.7) 
Introducing Fourier transforms in the x-space, (3.7) become 
Рат, (р, т) = 0, (42р? + 1 r) (p, A =A (3.77) 
and for a particle at rest- (p = 0, Pı= iM) 
Mry = 0, А ИО (9-795) 


hence it is seen that for p = 0 the field quantity is zero unless r, = 0 and r = 
= V72— 4? | M |. In this way the range of variability of г is restricted to polar 
angles Ф, 9 on a sphere with radius 22 | М |. 


4. Main features of the mass spectra 


We assume the viewpoint that the experimentally observed rest masses of 
particles consist of two parts of quite different origin. One of them (given by the 
eigenvalues of the mass operators appearing in the interaction free bilocal field 
equations) may be called the proper mass, the other (due to the interaction) may 
be called the field mass. It is to be expected that the mass of the electron is of 
a purely electromagnetic origin. From the ratio of the electromagnetic coupling 
constant to the nuclear forces coupling constant it is to be expected that the 
corresponding field masses due to nuclear coupling are larger than the electron 
mass roughly by two orders of magnitude. Thus, the investigation of the mass 
eigenvalues of the nucleon and pion families without taking account of their in- 
teraction will yield only a rough picture of the spectra. Conversly, taking account 
of the interaction will yield a fine structure of mass spectra with mass terms 
shifted up to 10? electron masses (m,). In this section we are concerned only with 
the main features of nucleon and pion mass spectra with neglect of the fine struc- 
ture caused by interaction. 

The nucleon family. The Fourier transform of the generalized Dirac equation 
written in terms of the Yukawa variables is 


PT 35 
(t aup mit) oon =0 (4.1) 
with ay, given by (2.12). For a particle at rest. (p = 0, p, = iM) we get 
a 9 
j = 2 
(„m + кд)» (M, r) = 0. (4 ) 


Going over to spherical coordinates in the r-space we have 


9f 9 9 3" Е 
EEE ae ‘= Ir! 4.3 
В, ar, = В, 5 +в (5 + : зи ( ) 
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with Bun; 

[= Вт b+ Ly (4.4) 


where т denotes the Pauli matrix-vector in the space of the f indices, / is the or- 
bital angular momentum in the r-space, в is a matrix commuting with / and anti- 
commuting with. Ва with the property =? = 1. в together with f, and, on the other 


hand, I may be represented as sets of independent two row and column matrices 


whose eigenvalues are + lfor e and (j + 5) for (т. 1+1) where j= 1,2, ad 


For the derivation of (4.3) see e.g. Dirac, Quantum Mechanics (1947). 
In consequence of the constraints introduced by the supplementary conditions 


д 1 
the terms with ——— and —— + — disappear, while г in (4.3) is to be replaced by 
Or, OT Paar 
the radius of the sphrere 72 | M |. Thus, 
ао! 
(= мж $ (+5 ) aM] yae, 0) - o (4.5) 


so that the mass eigenvalues M jare given by 


А ОЬ Ж наді tiger 
асту NRE р . 46) 


The lowest mass value should correspond to the nucleon which offers а possibility 
of determining the universal length or rather the ratio as follows 
A 


Va = 7,4. 10-15 cm. (4.7) 


The ratio of the two lowest mass terms 


M. 
i Do 
AE үз (4.8) 


is too large to interprete M,, as the mass of the Г, particle (which is about 1,2 
times larger than the nucleon mass). However, a better agreement could not be 
expected without taking account of the fine structure. 


The pion family. The second order wave equation satisfied by pions is 


9? 9? 
(а 3 — 4a? 2 a (x, г) =0 (4.9) 


witha, given by (2.12). Going over to Fourier transforms v (p. ) and denoting 
р? = —M?, we have 


92 / 
c =) y(p,r) = 0. | (4.10) 
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Introducing spherical coordinates and taking account of the constraints, the terms 
with derivatives with respect to r, and r disappear, and we are left with 


nca tn 


so that the mass eigenvalues are given by 


4 IST 
qu 2 == 
Матаны» алаш Шаа (4.12) 


The lowest state is the л meson with M, = 0. This is to be interpreted as follows: 
the z meson possesses no proper mass but only a field mass (which may well be 
of the order 10? m, on account of a large value of the nuclear forces coupling 


constant). The next terms. are 
май 
10 


4 
АА 


sete. (4.13) 
Жат À 


Assuming the value (4.7) previously determined for the nucleon, we get 
М, = 880 т,, М, = 1160 m,, etc. ... (4.14) 


corresponding roughly to the masses of t-meson and heavier mesons of the pion 
family. Better agreement cannot be expected without taking account ‘of inter- 
action. 

We notice a qualitative agreement with experiment: the formalism accounts 
for the fact that the lightest pion has rather a small mass while the next members of 
the family possess masses of the order of magnitude of 1000 m,. The situation is 
quite different with nucleons possessing masses of the order of 2000 m, and 
higher. 


5. Fine structure of the mass levels 


In consequence of interaction every particle acquires a field mass so that 
the wave equation (of a renormalizable field theory) may be written for one particle 
states simply as an interaction-free equation with an additional mass x. It should 
be stressed that x is not a new arbitrary constant introduced "at hoc", but is 
determined uniquely for each type of particles by the character of interactions . 
and by the magnitude of the coupling constants. In the framework of a non-local 
(bilocal) theory it is a well defined, finite constant. | | | 

The nucleon family. For the nucleon we get a generalized Dirac equation 


with field mass 


| 9 а д 
SES RE ns 3 — 3 =0 5.1 
б 9x, 2 get Pa Ir, D АЗ Sa 
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or, for a particle at rest, 


( a,M bs iB, —— 2-4) v(M,r) = 0. (5.17) 


Taking account of the haat ah we * the following equation for the mass 


eigenvalues 
ta j+} ‚el 
M + 4=0, = —,—.... 52 
+М + ИТ О ИККО (5.2) 


This yields four possible families of particles and antiparticles with absolute 
values of the masses 


2|МД- E Ere 2)). (5.2) 


8 а OR 
It is to be expected that x is much smaller than m so that a negative sign before a 


cannot be realized. Thus, in consequence of the existence of the field mass, each 
level j = 4, 2, ... splits into two (but not four) separate levels given by 


| M. sede х P 
ЗЕ TA (але? 1 + 5 (5.3) 


with equidistant splitting 


M,,—M,_ =x ; (5.4) 

We are unable yet to give a theoretical computation of the field mass but may 

derive its value from experiment. It is to be expected that the two lowest levels M АЕ 

апі М, „correspond to the nucleon (1840 m,) and to the V, particle (2185 m,) 
2 


respectively, whence 


x = 345 m, (5.5) 
This semi-empirical result enables us to compute a corrected value of = from 
MEL. 4i a Nr 1 1840 C 
Z 81? 4 yrs RETE t 
which yields 
a A je 
Eyre 33,79 and Ves 6,79 RIO tom | (5.7) 


The mass difference between the Г, particle and the nucleon being used to deter- 


mine the value of the field mass, the higher. mass levels may be computed by 
(5. 3) . We find 


oe 2670 m, 7 Мет 3015 m, , etc. ... (5.8) 


The value 2670 m is consistent with the recent findings of a charged particle A 
decaying into a neutral V, and a charged л. Our formalism yields for this process 
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a release of energy of Q = 108 MeV. This value is close to that observed by means 
of photographic plates (and referred to as a decay of the J particle), so that we 
are inclined to identify particle J with A. This interpretation requires the existence 


of a selection rule forbidding the transitions of the type 
= 


М, > M, (5.9) 


Most probably only transitions between a state + and a state — are allowed, 
so that the particle M, A would decay into anucleon (M, ) with emission ofa few 
5 $- 


pions (during a very short lifetime). 


The pion family. The second order wave equation corrected for the field 
mass is 


( 9? 4 9? 

Дал 2) v(x, г) = 0, (5.10 
9х2 Or, ) 
whence 

ii +1 
M? — 4а? ® — X? = 0) (5.11) 


with a? given by (2.12). The mass eigenvalues are 


Ee isa Тыл 
La 104 ЕО ОХ erm d om 5 == Л у (5.12) 


The lowest mass eigenvalue is identical with the field mass : 
M, = х =273m. і (5.13) 


In consequence of the above semi-empirical determination of the field mass, 
we get corrected values of the higher order masses from the formula 


M, 4 a? 3 1 |р 
EAE "we AU Ut ner. 5 14 
х ТЕ. Даҳ IOs brin 79 5 (6.14) 
Denoting the field masses (5.5) and (5.13) by xy and x,, we find from (5.7) 

2 2 4 
eee ( =) = 186 300, (5.15) 
даха А ж„ 
whence 

M, = 3,566 Му ог М, -974т,. (5.16) 


In the same way we obtain | 

М, = 1270 т, М» = 1506 m,, etc. ... (5.17) 
The above values agree, within experimental errors, with т,х and x mesons res- 
pectively. The 4 mesons (~550 m,) are not found among the members of this 
family. It may be that they are not elementary particles but merely bound states 


of two л mesons.  ; | 
I would like to thank Dr Yukawa and Dr Pais for letting me see their papers 


before publication. 


as" c | J. Rayski 


KPATKOE COJIEPZKAHME 


Й. Раиски, Спектр масс элементарных частиц. 

Доказано, что обобщенные уравнения поля (уравнения Шредингера, Гор- | 
дона и Дирака) с придаточным условием взаимности, требованным в рамках 
нелокальной (билокальной) теории поля, ведет к сквантованию массы. Имеем 
возможность классифицировать семьи частиц согласно (полной или MOMO- 
винчатой) величине спина и спина изготопического. Поданы спектры масс è 
семьи ядерных частиц и семьи мезона л Массы состоят из двух частей; 
обыкновенной массы и массы полевой. Появление полевой массы ответственно 
за тонкое строение уровней масс. : | 
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